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ABSTRACT OF THE DISSERTATION 



A Search for New Physics with High Mass Tau Pairs in 
Proton-Antiproton CoUisions at = 1.96 TeV at CDF 

by Zongru Wan 
Dissertation Director: Professor John Conway 

We present the results of a search for new particles decaying to tau pairs using the data 
corresponding to an integrated luminosity of 195 pb~^ collected from March 2002 to 
September 2003 with the CDF detector at the Tevatron. Hypothetical particles, such 
as Z' and MSSM Higgs bosons can potentially produce the tau pair final state. We 
discuss the method of tau identification, and show the signal acceptance versus new 
particle mass. The low-mass region, dominated by Z — rr, is used as a control region. 
In the high-mass region, we expect 2.8 it 0.5 events from known background sources, 
and observe 4 events in the data sample. Thus no significant excess is observed, and we 
set upper limits on the cross section times branching ratio as a function of the masses 
of heavy scalar and vector particles. 
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Chapter 1 
Introduction 

The Standard Model (SM) combines the electroweak theory together with Quantum 
Chromodynamics (QCD) of strong interactions and shows good agreement with colhder 
experiments. However the SM does not include gravity and is expected to be an effective 
low-energy theory. 

The Fermilab Tevatron is currently the high energy frontier of particle physics and 
delivers proton-antiproton collisions at high luminosity. 

The Run II of the Collider Dectector at Fermilab (CDF) continues the precision 
measurements of hadron collider physics and the search for new physics at and above 
the electroweak scale. With the precision capability at the energy frontier, we can 
attack the open questions of high energy physics from many complementary directions, 
including: the properties of top quark, the precision electroweak measurements, e.g. 
mass of the W boson, the direct searches for new phenomena, the tests of perturbative 
QCD at Next-to-Leading-Order and large Q^, and the constraint of the CKM matrix 
with high statistics of the B decays. 

This thesis is about a direct search for new particles decaying to tau pairs. The 
evidence for such new particles is that at accessible energies the events with tau pairs 
deviate clearly and significantly from the SM prediction. 

In Run I CDF recorded an unusual event in which there were two very high energy 
T ^ hu candidates nearly back-to-back in direction. Figure im shows a display of the 
event. This event was recorded in the data sample from the missing transverse energy 
trigger, and was noticed in the context of the Run I charged Higgs search pp. In Run I, 
a posteriori, it was not possible to estimate a probability for observing such an event, 
though less than about 0.1 such events were expected from backgrounds, including 
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E=l 60 GeV 




Figure 1.1: Run I high-mass di-tau candidate event. The left plot shows energy mea- 
surement in calorimeters and the event is very clean. The right plot shows the display 
in the transverse plane. The three-prong identified tau object has energy 160 GeV. The 
one-prong identified tau object has energy at least 135 GeV. There is also a significant 
missing transverse energy indicating significant neutrinos. The scale of the invariant 
mass of the two tau objects and the neutrinos is above 300 GeV/c^. 

^ rr Drell-Yan {qq Z/7* ^ /+/~). 

Various new physics processes can lead to very high-mass tau pairs, for example, the 
new vector boson Z' tt predicted in the extension to the Standard Model by adding 
a new U(l) gauge symmetry and the pseudoscalar Higgs boson A ^ tt predicted in 
the minimum supersymmetric extension of the Standard Model (MSSM). The known 
backgrounds are from the high-mass tail of Drell-Yan processes (mainly Z/^* tt) as 
well as jet— > r fakes from VF-l-jets, QCD di-jet, and mutli-jet events. 

In this analysis we search for such signal processes by performing a counting exper- 
iment. We select events with e + T^, fi + T^, and T^ + Th (here, "r/i" means a r hadronic 
decay). We construct an invariant mass which we call rrims using the four- vector sum 
of the lepton, the tau, and the missing transverse energy vector (ignoring in the latter 
the z component). The region which has m^j^ > 120 GeV/c^ is defined as the signal 
region, while the region which has rriyis < 120 GeV/c^ is retained as a control region. 
We perform a blind analysis in the signal region, i.e., we do not look at the data in the 
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signal region until we have precisely estimated the backgrounds. If there is a significant 
excess over the known backgrounds, we have discovered new physics; otherwise, we set 
limits on the possible signal rates. 

The thesis is organized as follows: theorectical models including the SM, extensions 
of the SM, and high-mass tau pair phenomenology are described in Chapter |21 The 
experimental appratus including the Fermilab Accelerator and CDF detector is intro- 
duced in Chapter |21 We discuss the logic behind the analysis in Chapter 0] Particle 
identifications for tau, electron and muon, and the study of missing transverse energy 
are discussed in detail in Chapter |SJ The data samples and event selection are dis- 
cussed in Chapter El The low-mass control region background estimate, uncertainties, 
and the observed events are discussed in Chapter[7| The high-mass signal region, signal 
acceptance, background estimate, and uncertainties are discussed in Chapter |H1 The 
results of the observed events after opening the box, and the method to extract limit 
are discussed in Chapter IHI Finally, the conclusion is presented in Chapter 1101 
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Chapter 2 
Theoretical Model 

The goal of elementary particle physics is to answer the following fundamental questions: 

• What is the world made of? 

• How do the parts interact? 

The Standard Model (SM) P of particle physics is a beautiful theory which attempts 
to find the simplest model that quantitatively answer these questions. The thousands 
of cross sections and decay widths listed in the Particle Data Group (PDG) [3^, and all 
of the data from collider experiments, are calculable and explained in the framework of 
the SM, which is the bedrock of our understanding of Nature. 

Building on the success of the SM, ambitious attempts have been made to extend it. 
This thesis is concerned about a direct search for new particles decaying to two taus. 
The phenomenology of tau pairs, namely the production rates of intermediate bosons 
and the branching ratio of their decays to tau pairs, in the framework of the SM and 
some of the extensions will be presented in this chapter. 

2.1 The Standard Model 

The SM elementary particles include the fermion matter particles and the force carriers. 
There are three generations of fermion matter particles: leptons and quarks. The second 
and third generations have the same quantum numbers of the first generation, but with 
heavier masses. The masses of the leptons and quarks are listed in Table 12.11 The 
force carriers include the gluon for the strong interaction, and the photon, the W and Z 
vector bosons for the electroweak interaction. The masses of the force carriers are listed 
in Table The Higgs boson predicted in the SM is a fundamental scalar particle and 
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Generation 


Particle 


Mass [GeV/c^] 


I 


electron neutrino 









electron 


e 


0.00051 




up quark 


u 


0.002 to 0.004 




down quark 


d 


0.004 to 0.008 


II 


muon neutrino 









muon 




0.106 




charm quark 


c 


1.15 to 1.35 




strange quark 


s 


0.08 to 0.13 


III 


tau neutrino 


Vr 







tau 


T 


1.777 




top quark 


t 


174.3 ± 5.1 




bottom quark 


b 


4.1 to 4.4 



Table 2.1: Three generations of leptons and quarks in the Standard Model and their 
masses. 



Force 


Carrier 


Mass [GeV/c^] 


electromagnetic 
charged weak 
neutral weak 
strong 


photon 7 
W boson 
Z boson 
gluon g 




80.425 ± 0.038 
91.1876 ± 0.0021 




Table 2.2: Force carriers in the Standard Model and their masses. 



has special interaction strength proportional to the mass of the elementary particles. 
Since it is not discovered yet, it is not listed in Table 

The SU(3)cxSU(2)2,xU(l)y structure of the leptons and quarks is shown in Fig. 12.11 
The quarks are arranged in triplets with respect to the color gauge group SU(3)c, with 
indices as red (r), green (g), and blue (6). 



Qr 



(2.1) 



The left- and right-handed fermions have different transformation properties under the 
weak isospin group S\J{2)l. The left-handed fermions are arranged in doublets, and 
the right-handed fermions are arranged in singlets. There is no right-handed neutrino 
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Figure 2.1: SU(3)cxSU(2)ixU(l)y gauge symmetries of fermions in the Standard 
Model. Quarks have three color degrees-of- freedom, while leptons are colorless. Left- 
handed fermions are arranged in SU(2) weak isospin doublets and right-handed fermions 
are arranged in SU(2) singlets. Each fermion also has U(l) weak hyper-charge. The 
interactions are uniquely specified by the gauge symmetries. 



in the SM. 

Leptons: 

Quarks: 




^J'R 



TR 



L (2.2) 

UR dR CR SR tR hR 



Table lT^ lists the transformation properties, i.e., the quantum numbers, of the fermions 
of the first generation under the gauge groups. The hypercharge of U(l)y is related to 
the electric charge by Q = T| + The assignments of the quantum numbers to the 
second and third generations are the same. A brief review about how this structure 
emerges is given in Appendix 1X1 

The interactions are uniquely specified by the SU(3)cxSU(2)/,xU(l)y gauge sym- 
metries. All of the gauge bosons and fermions acquire mass by the Higgs mechanism 
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Q 


r| 


Y 


C 







1/2 


-1 







-1 


-1/2 


-1 





eR 


-1 





-2 





UL 


2/3 


1/2 


1/3 


r,g,b 


dL 


-1/3 


-1/2 


1/3 


r,g,b 


UR 


2/3 





4/3 


r,g,b 


dR 


-1/3 





-2/3 


r,g,b 



Table 2.3: Quantum numbers of the fermions. 





Left Coupling 


Right Coupling 


Higgs H ^ ff 


V 


'inj_ 

V 


Strong g ^ qq 




93 \a 


EM 7 ^ // 


eQf 


eQf 


Weak ff 
lui 
qq' 


92 
V2 

V 


-^^^^'GwQf 





Table 2.4: Couplings to fermions in the Standard Model. 



It introduces an extra Higgs boson, and its physical vacuum is spontaneously broken in 
the field space of the Higgs potential. The quark states in charged weak interactions me- 
diated by bosons are not the physical states, but rather a quantum superposition of 
the physical states, described by the CKM (Cabbibo-Kobayashi-Maskawa) matrix 



V / weak 



Vud Vus Vub 
Vcd yes Vcb 

Vtd Vts Vtb 



(2.3) 



mass 



The topic of this thesis is mostly related to the fermion couplings. The couplings 
to fermions in the SM are listed in Table 12.41 A very detailed review with explicit 
derivations on these topics starting from the gauge symmetry to the couplings to the 
fermions in the SM is given in Appendix lEl 

In spite of its tremendous success in explaining collider results, there are still many 
unexplained aspects in the SM. The set of group representations and hypercharge it 
requires are quite bizarre, and there are 18 free parameters which must be input from 
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experiment: 3 gauge couplings (usually traded as e, sin^ 9w and g^), 2 Higgs potential 
couplings (usually traded as niz and mn), 9 fermion masses, and 4 CKM mixing pa- 
rameters. Do particle masses really originate from a Higgs field? Can all the particle 
interactions be unified in a simple gauge group? What is the origin of the CKM ma- 
trix? The ultimate "theory of everything" should explain all of these parameters. The 
imaginary goal, for example, is probably to express everything in terms of the Planck 
constant h, the speed of light c, the mathematical constant vr, and without any free 
parameters. That would be an amazing accomplishment. There are still many things 
to do in particle physics in the direction to find the simplest model and many exciting 
challenges are ahead! 

2.2 Extensions to the Standard Model 

One interesting extension to the SM is to add a new U(l) gauge group. This predicts 
a new Z' gauge boson |S] at high energy scale. We will use the Z' as our model to 
calculate the signal acceptance for any kind of new vector boson. 

Another interesting extension is supersymmetry ,7 , which is motivated by the desire 
to unify fermions and bosons, shown in Fig. 12.21 For each fermion (lepton and quark) 
it predicts a bosonic super partner (slepton and squark), and for each gauge boson 
it predicts a fermionic super partner (gaugino). There is a divergence from scalar 
contributions to radiative corrections for the Higgs mass in the SM, while the new 
fermion loops appearing in supersymmetry have a negative sign relative to the scalar 
contributions, thus cancel the divergence. We will use the pseudoscalar Higgs particle 
A, one of the Higgs particles predicted in the minimal supersymmetric extension of the 
Standard Model (MSSM) !8_ as our model to calculate the signal acceptance for any 
kind of new scalar boson. 

2.3 High Mass Tau Pairs 

At the Tevatron, the tau pair production in the SM is through the Drell-Yan process, 
PP 7* — > TT, as shown in Fig. 12.31 The center-of-mass energy oipp collisions at the 
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Quark 




Slepton 

Figure 2.2: Particles in the Supersymmetry Theory. 

Tevatron is 1.96 TeV. At the parton level, one incoming quark from a proton and the 
other anti-quark from an anti-proton collide via an intermediate boson which decays 
to two outgoing taus. The details about how to calculate cross sections are shown in 
Appendix [0 and the mass spectrum of the final two taus is shown in Fig. 12.41 We 
perform a direct search for new hypothetical particle in high mass region by its decay 
to two taus X — > TT. The low mass region of the SM processes ^* /Z — > rr is the 
control region and its high mass Drell-Yan tail is the major background for this search. 

The two extensions described above are shown in Fig. 12.51 For U(l) extension, 
we consider the simplest model with the same interactions as the Z boson in the SM, 
called the sequential Z', and the only unknown parameter is the mass of the new 
gauge boson. The MSSM requires two Higgs doublets and the ratio of the two Higgs 
expectation values is defined as tan /?, which is undetermined and should be treated as 
a free parameter. Thus the A boson is governed by one more free parameter in addition 
to its mass. 

The couplings to fermions in the SM are listed in Table ITU For each mass point 



10 




Figure 2.3: Tau pair production pp ^ 7*/Z — > rr in the Standard Model. 




100 200 300 400 500 600 
m,, [GeV/c^] 

Figure 2.4: High mass tau pair search. Low-mass region including the Z peak is the 
control region. High-mass region is the signal region. The high-mass tail of the Drell- 
Yan process is the main background of this search. The signature of new particles is a 
sigficant deviation from the known backgrounds, such as X shown in this plot. 
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of the sequential Z', we can use the same couphngs to fermions as the Z boson in the 
SM and repeat the procedure to calculate the cross section. The leading order cross 
section ctq is subject to a correction K factor 9_ such that the corrected cross section 
(T = (1 + correction) x uo = ivT x a^. Including the K factor, the predicted cross section 
versus mass for the sequential Z' is shown in Fig. 12.61 

The SM requires one Higgs doublet with a coupling of the SM Higgs boson to 
fermions as m//u, where m/ is the fermion mass and v is the vacuum expectation value 
of the SM Higgs boson, about 246 GeV. Therefore Higgs boson prefers to couple to 
the fermions in the heaviest generation. In the MSSM, at large tan/3, the coupling 
of ^ —> rr and A ^ hb are enhanced to mftanP/v, whereas the coupling of ^ — > 
is suppressed to nit cot P/v when the top quark is kinematically available, i.e. niA > 
2mt ^ 350 GeV/c^. We use the programs HIGLU [TOI and HDECAY [1^ to calculate 
the next-to-leading-order cross section oi gg — > A — > rr. They are also shown in 
Fig.EEl 
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Figure 2.5: Tree-level Feynman diagrams for the productions at pp collider and decays 
of Z' predicted in U(l) extension and pseudoscalar A predicted in minimum supersym- 
metric extension of the Standard Model. 




Figure 2.6: Theoretical signal a{pp — > X) ■ B{X tt). 
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Chapter 3 

The Tevatron Accelerator and the CDF Detector 

Fermilab is the home of the highest energy particle accelerator in the world, the Teva- 
tron. The center-of-mass energy of proton- antiproton (pp) collision is y/s = 1.96 TeV. 
We shall describe the Tevatron accelerator and the Collider Detector at Fermilab (CDF) 
in this chapter. 

3.1 Fermilab's Accelerator Chain 

Protons and antiprotons have equal and opposite electric charge. The advantage of 
pp collider is that p and p travel in opposite directions through the magnets and a pp 
collider can be built with one ring of magnets instead of two. The disadvantage is that 
it is difficult to produce and accumulate p at a high efficiency. 

The aerial view of Fermilab is shown in Fig. 13.11 The Fermilab's accelerator chain 
is shown in Fig. 13.21 It consists of the Proton/ Antiproton Sources (8 GeV), the Main 
Injector (150 GeV), the Recycler, and the Tevatron (980 GeV). 

The Proton Source includes the Cockcroft- Walton, the Linear Accelerator (Linac), 
and the Booster. The Cockcroft- Walton uses DC power to accelerate H~ ions to 750 
KeV. The Linac uses Radio Frequency (RF) power to accelerate H~ ions to 400 MeV. 
The electrons are stripped off and the bare protons are injected into the Booster. The 
Booster uses RF cavities to accelerate protons to 8 GeV. 

The Anti-proton Source includes the Target Station, the Debuncher and the Accu- 
mulator. A bunched beam of 120 GeV protons from the Main Injector hits a Nickel 
Target to make anti-protons and other particles as well. The particles are focused 
with a lithium lens and filtered through a pulsed magnet acting as a charge-mass spec- 
trometer to select anti-protons. The antiproton beam is bunched since the beam from 
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Figure 3.1: Aerial view of Fermilab showing the Main Injector in the foreground, the 
Tevatron colhder ring and the fixed target facihties in the background. 




Figure 3.2: Fermilab's accelerator chain consists of the 8 GeV proton source, the 8 
GeV anti-proton source, the Main Injector, the Recycler for recycling the precious anti- 
protons, and the Tevatron. The Main Injector accelerates protons and anti-protons to 
150 GeV. The Tevatron ramps up their energies to 980 GeV. The center-of-mass energy 
of pp collision is thus 1.96 TeV. The linear accelerators for the fixed targed experiments 
are also shown. 
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the Main Injector is bunched and the antiprotons have a wide range of energies, po- 
sitions and angles. The transverse spread of the beam out of the Target Station is 
"hot" , in terms analogous to temperature. Both RF and stochastic cooling systems are 
used in the momentum stacking process. The Debuncher exchanges the large energy 
spread and narrow time spread into a narrow energy spread and large time spread. The 
Accumulator stacks successive pulses of antiprotons from the Debuncher over several 
hours or days. For every million protons that hit the target, only about twenty 8 GeV 
anti-protons finally get stacked into the Accumulator. 

Protons at 8 GeV from the Booster are injected into the Main Injector. They are 
accelerated to 120 GeV for fixed target experiments or 150 GeV for injection into the 
Tevatron. Antiprotons at 8 GeV from either the Accumulator or the Recycler are 
accelerated to 150 GeV in the Main Injector and then injected into the Tevatron. 

The Recycler is placed directly above the Main Injector beamline, near the ceiling. 
One role of the Recycler is a post-Accumulator ring. Another role, and by far the 
leading factor in the luminosity increase, is to act as a recycler for the precious antipro- 
tons left over at the end of Tevatron stores. It is a ring of steel cases holding bricks 
of "refrigerator" magnets (the same permanent magnet used in home refrigerators). 
Permanent magnets do not need power supplies, cooling water systems, or electrical 
safety systems. The Recycler is a highly reliable storage ring for antiprotons. 

The Tevatron was the world's first superconducting synchrotron. A magnet with su- 
perconducting coils has no electrical resistance, and consumes minimal electrical power, 
except that is needed to keep the magnets cold. The particles of a beam are guided 
around the closed path by dipole magnetic field. The radius of the circle is 1000 meters. 
As the beam energy is ramped up by RF cavities from 150 GeV to 980 GeV, the bend- 
ing magnetic field and the RF frequency must be synchronized to keep the particles in 
the ring and this enables a stable longitudinal motion. The stability of the transverse 
motion is achieved with a series quadrupole magnets with alternating gradient. 

Luminosity is a measure of the chance that a proton will collide with an antiproton. 
To achieve high luminosity we place as many particles as possible into as small a collision 
region as possible. At the interaction point, the two beams of p and p are brought 
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Figure 3.3: pp collision. 



together by special quadrupole magnets called Low Beta magnets, shown in Fig. 13.31 
The current status (at the writing of the thesis) of the luminosity is shown in Fig. 13. 4| 
and the integrated luminosity delivered and to tape is shown in Fig. 13.51 

The design value for the peak instantaneous luminosity during Run II is 2 x 10^^ 
cm^^s^^. Typically a year allows 10^ seconds of running at the peak instantaneous 
luminosity. This is about one third of the actual number of seconds in a year, which 
accounts both for the drop in luminosity and for a normal amount of down-time. Using 
the conversion constant 1 fb = 10~^^ cm^, the design value corresponds to an integrated 
luminosity about 2 fb~^ per year. Ultimately it is hoped that an integrated luminosity 
of 8—10 fb^^ can be attained in Run II. The total number of events in a scattering 
process is proportional to the luminosity and the cross section a of the process, 

N = La (3.1) 

We can get a rough sense of the reach for new physics and the challenge of enhancing 
signal and suppressing background by considering the following examples. At a center- 
of-mass energy of 1.96 TeV, we have 

cr(pp anything) ~ 75 mb (3.2) 
a{pp — > + anything) 6 pb (3-3) 
a(j)p — > hZ + anything) 75 fb (3-4) 
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Figure 3.4: Run II instantaneous initial luminosity. 
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Figure 3.5: Run II integrated luminosity. 
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3.2 The CDF Dectector 

The CDF detector is cylindrically symmetric around the beamhne. A sohd cutaway 
view is shown in Fig. 13.61 and an elevation view is shown in Fig. 13.71 It is a general- 
purpose solenoidal detector with tracking system, calorimetry and muon detecion. The 
tracking system is contained in a superconducting solenoid, 1.5 m in radius and 4.8 m 
in length. The magnetic field is 1.4 T, parallel to the beamline. The calorimetry and 
muon system are outside the solenoid. These sub-systems will be described in more 
details below. 

3.2.1 CDF Coordinate System 

The origin of the CDF detector is its geometric center. The luminous region of the beam 
at the interaction point has Gaussian profiles with ((Ti:, Uy, cJz ) beam ~ (0-003, 0.003, 30) 
cm. The pp collision point is not necessarily at the origin. 

The CDF detector uses a right-handed coordinate system. The horizontal direction 
pointing out of the ring of the Tevatron is the positive x-axis. The vertical direction 
pointing upwards is the positive y-axis. The proton beam direction pointing to the east 
is the positive z-axis. 

A spherical coordinate system is also used. The radius r is measured from the center 
of the beamline. The polar angle is taken from the positive z-axis. The azimuthal 
angle (p is taken anti-clockwise from the positive x-axis. 

At a pp collider, the production of any process starts from a parton-parton interac- 
tion which has an unknown boost along the z-axis, but no significant momentum in the 
plane perpendicular to the z-axis, i.e. the transverse plane. This makes the transverse 
plane an important plane in pp collision. Momentum conservation requires the vector 
sum of the transverse energy and momentum of all of the final particles to be zero. The 
transverse energy Et and transverse momentum pT are defined by 



Ex 
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Figure 3.6: Solid cutaway view of CDF II detector. 




Figure 3.7: Elevation view of CDF II detector. 
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Hard pp head-on collisions produce significant momentum in the transverse plane. 
The CDF detector has been optimized to measure these events. On the other hand, 
the soft collisions such as elastic or diffractive interactions or minimum-bias events, and 
by-products from the spectator quarks from hard collisions, have most of their energy 
directed along the beampipe, and will not be measured by the detector. 

Pseudorapidity rj is used by high energy physicists and is defined as 



Consider occupancy in a sample of large amount of pp collision events. Typically, 
particles in a pp collision event tend to be more in the forward and backward regions 
than in the central region because there is usually a boost along the z-axis, which could 
be shown in 6 occupancy of the particles of the events in the sample. Now we transform 
io rj. The derivative of r/ is 



A constant rj slice corresponds to variant 9 slice which is smaller in the forward and 
backward regions than in the central region. This can make the r] occupancy more 
uniform than 9 occupancy. For example, calorimeters are constructed in rj slices, instead 
of 9 slices. 

3.2.2 Tracking 

The tracking volume is surrounded by the solenoid magnet and the endplug calorime- 
ters as shown in Fig. K-{.81 The tracking system records the paths of charged particles 
produced in the pp collisions. It consists of a silicon microstrip system with radius 
from r = 1.5 to 28 cm and \rj\ < 2, and an open-cell wire drift chamber called central 
outer tracker (COT) j2] with radius from r = 40 to 137 cm and \rj\ < 1. 

The silicon microstrip is made from Si with a p-n junction. When p-type semi- 
conductors and n-type semiconductors are brought together to form a p-n junction, 
migration of holes and electrons leaves a region of net charge of opposite sign on each 
side, called the depletion region (depleted of free charge carriers). The p-n junction can 
be made at the surface of a silicon wafer with the bulk being n-type (or the opposite 
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LAYER 00 SVX II INTERMEDIATE SILICON LAYERS 

Figure 3.8: CDF II tracking volume. 

way) . By applying a reverse-bias voltage we can increase the depletion region to the full 
volume of the device. A charged particle moves through this depletion region, creates 
electron-hole pairs which drift and are collected at the surfaces. This induces a signal 
on metal strips deposited on the surface, connected to readout amplifiers. 

The silicon microstrip detector consists of three components: the Layer 00, the 
Silicon VerteX detector II (SVX II), and the Intermediate Silicon Layers (ISL). An 
end view is shown in Fig. 13.91 Layer 00 is physically mounted on and supported by 
the beam pipe. The sensors are single-sided p-in-n silicon and have a pitch of 25 /xm. 
The next five layers compose the SVX II and are double-sided detectors. The axial 
side of each layer is used for r-cp measurements and the sensors have a strip pitch of 
about 60 jjivn. The stereo side of each layer is used for r-z measurements. Both 90° 
and small-angle stereo sensors are used in the pattern (90, 90, —1.2, 90, -1-1.2) degrees 
and have a strip pitch of (141, 125.5, 60, 141, 60) /^m from the innermost to outermost 
layers. The two outer layers compose the ISL and are double-sided detectors with a 
strip pitch of 112 /im on both the axial and the 1.2° stereo sides. This entire system 
allows charged particle track reconstruction in three dimensions. The impact parameter 
resolution of SVX II -|- ISL is 40 /zm including 30 /im contribution from the beamline. 
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Figure 3.9: Silicon system. Figure 3.10: COT superlayers. 

The zq resolution of SVX II + ISL is 70 //m. 

The COT is arranged in 8 superlayers shown in Fig. 13.101 The superlayers are 
alternately axial and ±2° stereo, four axial layers for r-(j) measurement and four stereo 
layers for r-z measurement. Within each superlayer are cells which are tilted about 
30° to the radial direction to compensate for the Lorentz angle of the drifting charged 
particles due to the solenoid magnet field. Each cell consists of 12 layers of sense wires, 
thus total 8x12 = 96 measurements per track. 

The COT is filled with a mixture of argone:ethane = 50:50 which determines the 
drift velocity v. A charged particle enters gas, ionizes gas and produces electrons. There 
is an electric field around each sense wire. In the low electric field region, the ionization 
electrons drift toward the sense wire. In the high electric field region within a few 
radii of the sense wire, there is an avalanche multiplication of charges by electron-atom 
collision. A signal is induced via the motion of electrons. By measuring the drift time 
t (the arrival time of "first" electrons) at sense wire relative to collision time to, we can 
calculate the distance of the hit D = vAt. 

A track is formed from a series of hits, fit to a helix. We can measure the curvature 
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Figure 3.11: Calorimeter tower segmentation in rj — (f) space. 



of a track C = 1/R and then calculate transverse momentum pT = 0.3RB, with pT, 
R and B in the units GeV/c, m, and T, respectively. The hit position resolution is 
approximately 140 ^m and the momentum resolution cr{pj')/p'x = 0.0015 (GeV/c)^^. 

3.2.3 Calorimetry 

The CDF electromagnetic and hadronic sampling calorimeters surround the tracking 
system and measure the energy flow of interacting particles up to |?7| < 3.64. They are 
segmented in rj and (j) with a projective "tower" geometry, shown in Fig. 13.111 

The energy measurement is done by sampling calorimeters which are absorber and 
sampling scintillator sandwich with phototude readout. When interacting with the ab- 
sorber, electrons lose energy by ionization and bremsstrahlung, and photons lose energy 
by the photoelectric effect, Compton scattering and pair production. Both electrons 
and photons develop electromagnetic shower cascades. The size of the longitudinal 
shower cascade grows only logarithmically with energy. A very useful cascade param- 
eter is the radiation length Xq which is the mean distance for the to lose all but 
1/e of its energy. For example, for a 10 GeV electron in lead glass, the maximum 
electromagnetic shower is at about QXq and the 95% containment depth is at about 
IGXq. Hadrons lose energy by nuclear interaction cascades which can have charged 
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pions, protons, kaons, neutrons, neutral pions, neutrinos, soft photons, muons, etc. It 
is much more comphcated than an electromagnetic cascade and thus results in a large 
fluctuation in energy measurement. In analogy to Xq, a hadronic interaction length A 
can be defined. Hadronic showers are much longer than the electromagnetic ones. 

The central calorimeters consist of the central electromagnetic calorimeter (CEM) |15j . 
the central hadronic calorimeter (CHA) 16 , and the end wall hadronic calorimeter 
(WHA). At approximately 6Xq in depth in the CEM, at which electromagnetic show- 
ers typically reach the maximum in their shower profile, is the central shower maximum 
detector (CES). The CEM and CHA are constructed in wedges which span 15° in az- 
imuth and extend about 250 cm in the positive and negative z direction, shown in 
Fig. 13.121 There are thus 24 wedges on both the +z and —z sides of the detector, for 
a total of 48. A wedge contains ten towers, each of which covers a range 0.11 in pseu- 
dorapidity. Thus each tower subtends 0.11 x 15° in r] x (p. CEM covers < |r/| < 1.1, 
CHA covers < \r]\ < 0.9, and WHA covers 0.7 < |r/| < 1.3. 

The CEM uses lead sheets interspersed with polysterene scintillator as the active 
medium and employs phototube readout, approximately WXq in depth, and has an 
energy resolution 13.5%/ \/Et © 2%, where © denotes addition in quadrature. The 
CES uses proportional strip and wire counters in a fine-grained array, as shown in 
Fig. 13.131 to provide precise position (about 2 mm resolution) and shape information 
for electromagnetic cascades. The CHA and WHA use steel absorber interspersed with 
acrylic scintillator as the active medium. They are approximately 4.5A in depth, and 
have an energy resolution of 75%/ ^/Et © 3%. 

The plug calorimeters consist of the plug electromagnetic calorimeter (PEM) jl7j . 
and the plug hadronic calorimeter (PHA). At approximately 6Xq in depth in PEM is 
the plug shower maximum detector (PES). Fig. 13.141 shows the layout of the detector 
and coverage in polar angle 36.8° > ^ > 3° (1.1 < \r]\ < 3.64). Each plug wedge 
spans 15° in azimuth, however in the range 36.8° > 6 > 13.8° (1.1 < \rj\ < 2.11) the 
segmentation in azimuth is doubled and each tower spans only 7.5°. 

The PEM is a lead-scintillator sampling calorimeter. It is approximately 21Xq in 
depth, and has an energy resolution of 16%/\/^ © 1%. The PES consists of two layers 
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Figure 3.12: CEM/CES/CHA wedge. 



Figure 3.13: CES strip and wire. 
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Figure 3.14: PEM/PES/PHA layout. 



Figure 3.15: PES U and V layers. 



of scintillating strips: U and V layers offset from the radial direction by +22.5° and 
—22.5° respectively, as shown in Fig. 13.151 The position resolution of the PES is about 
1 mm. The PHA is a steel-scintillator sampling calorimeter. It is approximately 7A in 
depth, and has an energy resolution of 74%/\/^ ® 4%. 



3.2.4 Muon Chambers 



The muon chambers are situated outside the calorimeters. In addition to the calorime- 
ters, the magnet return yoke and additional steel shielding are used to stop electrons. 
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Figure 3.16: Muon stub matching to a track. 

photons and hadrons from entering the muon chambers. The muon is a minimum ion- 
izing particle which loses very little energy in detector materials. The muon's lifetime 
is long enough to allow it to pass through all the detector components, reach the muon 
chambers, and decay outside. 

A muon chamber contains a stacked array of drift tubes and operates with a gas 
mixture of argon:ethane = 50:50. The basic drift principle is the same as that of the 
COT, but the COT is a multi-wire chamber, while at the center of a muon drift tube 
there is only a single sense wire. The sense wire is connected to a positive high voltage 
while the wall of the tube is connected to a negative high voltage to produce a roughly 
uniform time-to-distance relationship throughout the tube. The drift time of a single 
hit gives the distance to the sense wire, and the charge division at each end of a sense 
wire can in principle be used to measure the longitudinal coordinate along the sense 
wire. The hits in the muon chamber are linked together to form a short track segment 
called a muon stub. If a muon stub is matched to an extrapolated track, a muon is 
reconstructed. This is shown in Fig. 13.161 

There are four independent muon detectors: the central muon detector (CMU) jTH], 
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Figure 3.17: Muon coverage in ry and (j). 

the central muon upgrade (CMP), the central muon extension (CMX), and the inter- 
mediate muon detector (IMU) . The muon coverage in rj — (p space is shown in Fig. 13.171 

The CMU is behind the central hadronic calorimeter and has four layers of cylindri- 
cal drift chambers. The CMP is behind an additional 60 cm of shielding steel outside 
the magnet return yoke. It consists of a second set of four layers with a fixed length in 
z and forms a box around the central detector. Its psuedorapidity coverage thus varies 
with the azimuth. A layer of scintillation counters (the CSP) is installed on the outside 
surface of the CMP. The CMU and CMP each covers < 0.6. The maximum drift 
time of the CMU is longer than the pp bunch crossing separation. This can cause an 
ambiguity in the Level 1 trigger (described in the next section) about which bunch the 
muon belongs to. By requiring CMP confirmation, this ambiguity is resolved by the 
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CSP scintillators. 

The CMX has eight layers and covers 0.6 < \r]\ < 1.0. A layer of scintillation 
counters (the CSX) is installed on both the inside and the outside surfaces of the 
CMX. No additional steel was added for this detector because the large angle through 
the hadron calorimeter, magnet yoke, and steel of the detector end support structure 
provides more absorber material than in the central muon detectors. The azimuthal 
coverage of CMX has a 30° gap for the solenoid refrigerator. 

The IMU consists of barrel chambers (the BMU) and scintillation counters (the 
BSU), and covers the region 1.0 < |r/| < 1.5. 

3.3 Trigger and Data Acquisition System 

The trigger system has a three- level architecture: level 1 (LI), level 2 (L2), and level 3 
(L3). The data volume is reduced at each level which allows more refined filtering at 
subsequent levels with minimal deadtime. The trigger needs to be fast and accurate to 
record as many interesting events as possible, while rejecting uninteresting events. 

Each sub-detector generates primitives that we can "cut" on. The trigger system 
block diagram is shown in Fig. 13.181 The available trigger primitives at LI are 

• XFT tracks, with cp and pT provided by the eXtreme Fast Tracker using the hits 
in the axial layers of the COT, 

• electrons, based on XFT and HAD/EM which is the ratio of the hadronic energy 
and the electromagnetic energy of a calorimeter tower, 

• photons, based on HAD/EM ratio, 

• jets, based on EM+HAD, 

• muons, based on muon hits and XFT, and 

• missing Et and sum Et which are the negative of the vector sum and the scalar 
sum of the energies of all of the calorimeter towers, respectively. 

The available trigger primitives at L2 are 
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Figure 3.18: Trigger system block diagram. 



• SVT, the Silicon Vertex Tracker trigger based on the track impact parameter of 
displaced tracks, 

• jet clusters, 

• isolated clusters, and 

• EM ShowerMax which is the strip and wire clusters in the CES. 

There are two important factors for trigger design: the time between beam crossing 
and N, the average number of overlapping interactions in a given beam crossing. 

We can have many bunches in the Tevatron to enhance the luminosity. Since the 
radius of the ring is 1000 m, a proton (or an anti-proton) at a speed very close to the 
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Figure 3.19: Average number of interactions per crossing for various bunches, as a 
function of instantaneous luminosity. 



speed of light circulates the ring once every 20 /is. To accomodate 36 bunches, the 
maximum bunch separation allowed is about 600 ns, and the Run Ila configuration is 
396 ns. The bunch separation defines an overall time constant for signal integration, 
data acquisition and triggering. 

Another key design input is the average number of overlapping interactions N, which 
is shown as a function of luminosity and the number of bunches in Fig. 13.191 [19]. For 
example, with 36 bunches, N is about 1 at 3 x 31 cm~^s~^ and about 10 at 4 x 32 
cm~^s^^. The trigger with fast axial tracking at LI can handle the former environment, 
but cannot handle the latter environment because of the presence of too many fake 
tracks. To be able to handle 4 x 32 cm^'^s^^ we would need 108 bunches and even that 
seems not enough, thus we will also need to upgrade the trigger to include, for example, 
stereo tracking at LI to suppress fake tracks. 
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Figure 3.20: Data flow of "deadtimeless" trigger and data acquisition. 



The data flow in the trigger system is constrained by the processing time, i.e. how 
fast a decision can be made to clear events at each level and the tape writing speed for 
permanant storage at the end of the triggering process. The implementation needs a 
sufficient buffer while filtering because any overflow means deadtime. The "deadtime- 
less" design for 132 ns crossing is shown in Fig. I3.2UI 

The LI decision occurs at a fixed time about 5.5 ^s after beam collision. LI is a 
synchronous hardware trigger. To process one event every 132 ns, each detector element 
is pipelined to have local data buffering for 42 beam crossings. The LI accept rate is 
less than 50 KHz which is limited by the L2 processing time. 



32 



The L2 decision time is about 20 fis. L2 is a combination of hardware and software 
triggers and is asynchronous. If an event is accepted by LI, the front-end electronics 
moves the data to one of the four onboard L2 buffers. This is sufficient to process a 
LI 50 KHz accept rate and to average out the rate fluctuations. The L2 accept rate is 
about 300 Hz which is limited by the speed of the event-builder in L3. 

L3 is purely a software trigger consisting of the event builder running on a large PC 
farm. The event builder assembles event fragments from LI and L2 into complete events, 
and then the PC farm runs a version of the full offline reconstruction code. This means 
that fully reconstructed three-dimensional tracks are available to the trigger decision. 
The L3 accept rate is about 75 Hz which is limited by tape writing speed for permanent 
storage. 

Once an event passes L3 it is delivered to the data-logger sub-system which sends 
the event out to permanent storage for offline reprocess, and to online monitors which 
verify the entire detector and trigger systems are working properly. 

The data used in this analysis were collected from March 2002 to September 2003. 
It was for 396 ns with 36 bunches and for luminosity about 3 x 31 cm~^s~"^. This means 
that the trigger (designed for 132 ns) was sufficiently capable to handle the timing of 
bunch crossing with no need to worry about multiple interactions in this environment. 
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Chapter 4 
Search Strategy 

This chapter describes the overall logic of the high-mass tau tau search. There are three 
steps: 

1. Use W ^ TV events to cross check the r identification efficiency. 

2. Use Z ^ TT events to study the low-mass control region with rrims < 120 GeV/c^. 

3. Examine the high-mass signal region with m^s > 120 GeV/c^ for evidence of an 
excess signalling new physics. 

Tau Hadronic Decays 

The dominant decays of r's are into leptons or into either one or three charged 
hadrons, shown in Table ETTl The following short-hand notations for r and its decays 
are used, 

Te T ^ evv (4.1) 
r ^ livu (4.2) 
Th T —> hadrons v (4.3) 

The leptonic decays cannot be distinguished from prompt leptons. So tau identification 
requires a hadronic tau decay only, with a mass less than 

m(r) = 1.777 GeV/c^ (4.4) 

The net charge of the charged tracks is ±1. But we will not cut on charge because for 
very high energy taus there is an ambiguity of charge sign for very straight tracks. 
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Decay Mode 


Final Particles 


BR 


Leptonic 




17.8% 






17.4% 


Hadronic 1-prong 




11.1% 






25.4% 






9.2% 






1.1% 






0.7% 






0.5% 


Hadronic 3-prong 




9.5% 






4.4% 



Table 4.1: Tau dominant decay modes and branching ratios. 

The characteristic signature of hadronically decaying taus is the track multiplicity 
distribution with an excess in the 1- and 3-track bins. The excess, about 2:1 in these 
bins, is related to the tau hadronic branching ratios to one or three charged pions. 
Quark or gluon jets from QCD processes tend not to have such low charged track mul- 
tiplicity, but have a broader distribution peaking at higher multiplicities (3-5 charged 
tracks). Other final particles, namely photons, electrons, and muons have mainly 0, 
1, or 1 tracks, respectively, which are different from tau hadronic decays too. Seeing 
the tau's characteristic track multiplicity signature is a very important indication that 
backgrounds are under control. 

Since a ■ B{W tv) is about ten times larger than a ■ B{Z ^ tt) j20j we will use 
W ^ TV events to cross check the tau identification efficiency. 

Di-Tau Visible Mass 

There are six final states for tau pairs, shown in Table 14.21 TeTg and t^t^ modes 
cannot be distinguished from the prompt ee or the prompt ^/x, respectively. TgT^ mode 
has a special signature, but its branching ratio is small and its final particles tend to 
have low energy. For this analysis, we will look for three golden final states with at 
least one hadronic decay. 
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Final States 


BR 


TeTh 


22% 




22% 


Thru 


41% 




3% 




6% 




6% 



Table 4.2: Tau pair final states and their branching ratios. 

The high-mass tau pair search will be based on just counting the number of events 
with some specified set of cuts. It is desirable to measure for some variable a distri- 
bution which agrees with the Standard Model in some range, but deviates from it in 
another, thus giving a more convincing signal while also providing an estimate of the 
new particle's mass scale. 

There are at least two missing neutrinos in the golden final states, and therefore six 
unknown momentum components. With only two constraints from the two components 
of the missing transverse energy and the two constraints from two tau masses, there is 
at least a 2-fold ambiguity. It is not possible to reconstruct the tau pair invariant mass 
in general. 

The mass of the sum of the two tau's visible momentum and the missing transverse 
energy with its z-component set to zero is called the visible mass, 

m^is = m(r^\, + r^.^ + $t) (4.5) 

The invariant mass of the irreducible Z ^ tt background peaks at m{Z) k. 91 
GeV/c^. The visible mass distribution will be broadened and peak at somewhere less 
than 91 GeV/c^. We will study the sample with rrims < 120 GeV/c^ for Z ^ tt cross 
check. After all of the cuts, we want the control sample to be dominated hy Z ^ 
TT background, with jet background under control and other backgrounds negligible. 
A successful cross check between data and MC in the low-mass region will give us 
confidence to go further to the high-mass region. 
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Blind Analysis 

If a new particle with high mass exists and the statistics are sufficient, it will show 
up in the high-mass signal region. The strategy we choose is a blind analysis. The 
data sample with nivis > 120 GeV/c^ will be put aside until all selection criteria are 
fixed and all backgrounds are determined. The principle of a blind analysis is to avoid 
human bias. If the selection cuts are decided by the distributions of high mass region 
in the real data sample, there will be a strong bias and the probabilities calculated are 
meaningless. Given good understanding of backgrounds, there will be two possibilities 
after examining the data in the signal region. Either one will observe a number of 
events statistically consistent with the expected background rate, or there will be an 
excess signalling new physics. 
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Chapter 5 

Particle Identification and Missing Transverse Energy 

High energy pp collisions can produce a large number of particles. As illustrated in 
Fig. 15.11 the CDF detector with its tracking system, calorimeter and muon chambers 
can identify the following particles by the following patterns: 

• photon: cascade showering in electromagnetic calorimeter, but no associated 
charged tracks; 

• electron: a track, and cascade showering in electromagnetic calorimeter; 

• muon: a track, minimum ionization energy deposit in calorimeter, and hits in 
muon chambers; 

• jet: an object which cannot be identified as an isolated photon, or an isolated 
electron, or an isolated muon is identified as a jet; 

• missing transverse energy {^t)'- an imbalance of transverse energy in the whole 
calorimeter. 

The final particles and the are reconstructed by CDF II offline programs. 
5.1 Monte Carlo Simulation 

Often we need to predict the output in the detector including the final reconstructed 
particles and the J^t of a particular interesting process and compare with data. Usually 
the phase space of an event of the pp collision is too complicated to be calculated 
analytically. In this case Monte Carlo (MC) simulation is used. It has become a 
powerful tool used in many research areas including high energy physics. 
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muon chambers 

Figure 5.1: Patterns for identifying photon, electron, muon, charged hadron, and jet. 
Neutrino induces missing transverse energy. 

A well-known MC example is the Buffon's Needle. It involves dropping a needle 
on a lined sheet of paper and determining the probability of the needle crossing one of 
the lines on the page. The remarkable result is that the probability is directly related 
to the value of the mathematical vr. Suppose the length of the needle is one unit and 
the distance between the lines is also one unit. There are two variables, the angle 9 at 
which the needle falls and the distance D from the center of the needle to the closest 
line. 9 can vary from 0° to 180° and is measured against a line parallel to the lines on 
the paper. D can never be more than half the distance between the lines. The needle 
will hit the line if < ^ sin 9. How often does this occur? The probability "P is 2/7r by 
integrating over 9. With a computer, we can generate a large sample of random needle 
drops. The probability V can be simply taken as the number of hits divided by the 
number of drops, yielding vr = 2/V. 

Here we discuss the basic techniques of MC simulation. For a one-dimensional 
integral, we can choose n numbers Xj randomly with probability density uniform on 
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the interval from a to b, and for each Xi evaluate the function f{xi). The sum of these 
function values, divided by n, will converge to the expectation of the function /. 



The central limit theorem tells us that the sum of a large number of independent 
random variables is always normally distributed (i.e. a Gaussian distribution), no 
matter how the individual random variables are distributed. To understand this, we 
can test with uniformly distributed random variable xi, X2, x^, X4, (a) xi is a uniform 
distribution; (b) xi + X2 is a triangle distribution; (c) xi + X2 + X3 is already close 
to a Gaussian distribution; (d) xi + X2 + X3 + X4 is almost like the exact Gaussian 
distribution. Applying this theorem, we know the MC method is particularly useful as 
we can also calculate an error on the estimate by computing the standard deviation. 



where an = (/^ — fn )^^^ and = ^J2i=i f^ixi). The convergence for numerically 
evaluating the integral goes as l/\/n with the number of function evaluation, n. And 
obviously if the distribution f{x) is flatter, then the (t„ is smaller for the same number 
of events in a sample generated. If there is a peak in the distribution such as the 
distribution of a resonance production, it is better to transform that variable to some 
other variable with a flatter distribution in order to converge faster. 

The generalisation to multi-dimensional integrals J f{x,y,z, ...)dxdydz... is straight- 
forward. We can choose n numbers of grid {x, y, z, ...) randomly with probability density 
uniform on the multi-dimensional phase space, and for each grid evaluate the function 
f{x,y,z,...). The sum of these function values, divided by n, will converge to the ex- 
pectation of the function /. A nice feature is that it will always converge as l/\/n, even 
for very high dimensional integrals. This can make the performance of the MC method 
on multi-dimensional integrals very efficient. 

In high energy physics, an event occurs with a probability in the phase space of 
the kinematic variables. A MC simulation generates a large number of random events 
according to the probability described by a model. With a large sample, we can get 




(5.1) 



{f{x)) = fn±^ 



(5.2) 
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Figure 5.2: Illustration of tau isolation cone definitions. 



the predictions of the model by looking at the distributions of the kinematic variables 
and the derived variables, and the correlations among the variables. By confronting the 
predictions with real data, it is possible to tell if a model describes Nature correctly. 

For this analysis, we use PYTHIA 6.215 program j2J with CTEQ5L parton density 
functions (PDF's) [221 generate the large samples of the processes of pp collision, 
such as ^ 7*/Z ^ rr, pp ^ Z' ^ tt, pp ^ A ^ tt, and use TAUOLA 2.6 [2S1 
to simulate tau decays. We use GEANT 3 |24j to simulate the response to the final 
particles in the CDF II detector. 

5.2 Tau Identification 

Tau leptons decay predominantly into charged and neutral pions and suffer from large 
backgrounds from jet production. Hadronic tau decays appear in the detector as narrow 
isolated jets. The most powerful cut to suppress the jet background is in fact isolation, 
requiring no other tracks or vr'^s near the tau cone. To do this we define a signal cone 
and an isolation cone around the direction of the seed track (the track with the highest 
Pt) and then require that there is no track or vr*^ between the signal cone and the 
isolation cone. This is shown in Fig. 15.21 
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5.2.1 Cone Size Definition 

There are two useful cone size definitions. One is to construct a cone in AR defined 
below which has relativity invariance under a boost along the z-axis. The other is to 
construct a cone in three-dimensional separation angle, a, which has geometry invari- 
ance. Below we discuss why AR is chosen as cone size definition for jet identification 
and why a is chosen as cone size definition for hadronic tau identification. 

We start with the discussion on relativity invariance. For a particle under a boost 
(3 = v/c along the z-axis and 7 = (1 — its four-momentum {p^i Py, Pz, E) is 

transformed to 



10 
10 

7 /?7 

V /37 7 ; 



(5.3) 



Py 
Pz 

\E J 

The Px and py components in the transverse plane are not changed, while the Pz com- 
ponent and the energy are changed. Rapidity is defined by 



Px 

Py 

j{p, + pE) 

\ ^{Ppz + e) J 



2 ""e-p. 



(5.4) 



Using tanh f3 



M - 1 In IM 



1-/3 



it is easy to check that rapidity has a nice additive property 



under the boost along the z-axis. 



y ^ y + tanh ^ /3 



(5.5) 



For ultra-relativistic particle with p S> m, we have Pz/E pz/p = cos 9. Using 
cos9 = (1 — tan^ |)/(l-|-tan^ |), the rapidity is well approximated by pseudorapidity rj, 



?7 = — In tan - 



(5.6) 



Particles in a jet deposits energy in the calorimeter towers. For the traditional 
cone jet algorithm, we can call the tower with Et above a seed threshold as the seed 
(abbreviated as s), and the other towers with Et above a shoulder threshold as shoulders 
(abbreviated as h). To identify a jet, we can put the seed at the center and make a 
cone starting at a reconstructed interaction vertex point and around the seed to include 
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the shoulders. Since the transverse components of a particle's four-momentum are not 
changed under the unknown boost (3 of the parton-parton system along the z-axis, (f) is 
not changed. For an ultra-relativistic particle, rj is a good approximation of its rapidity. 
We have 

^ ^^^^ 

Vs 118 + tanh (3, rih rjh + tanh /? 

The separations in (p and rj are not changed under the unknown boost along the z-axis, 

A<p = cl)h-(ps ^ A<^ ^^^^ 

^V = Vh-Vs Ar/ 
Therefore the separation in AR which is constructed in the combination of Acj) and A77 
is not changed under the unknown boost along the z-axis, 

AR = ^{AriY + {A4>Y ^ AR (5.9) 

Given the Et and the configuration (shape) of a jet, whatever the magnitude of the 
boost along the z-axis of the parton-parton system is, or, equivalently, whatever the 
direction of the seed of the jet is, we can use the same cone to include or exclude a 
tower into the jet by calculating its separation in AR to the seed. Thus AR is a very 
useful shape variable for jet identification. 

It also makes sence that there is a strong correlation between the two variables Et 
and AR: a higher Et should give a smaller cone in AR to include all of the final 
particles, e.g. of a jet. It is very common that there arc hundreds of final particles after 
a pp collision. The problem is that the energy of a jet in real data cannot be measured 
before a cone is actually constructed, otherwise there is no constraint to tell which 
tower should be included or excluded. Jet identification usually starts with a large and 
constant cone around a seed. The towers with significant energy in the cluster may or 
may not be contiguous. The energy of the jet is determined afterwards by summing up 
the energies of all of the towers in the cluster. 

Now consider hadronic tau identification with a narrow cone and small number 
of final particles. The situation is quite different from jet identification. Since there 
are only a small amount of final particles, each final particle has significant energy. 
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And since all of the final particles are in a narrow cone, they make a narrow and 
contiguous cluster with significant energy in each tower. This constraint of a narrow and 
contiguous cluster with significant energy in each tower tells us that we can determine 
energy first, and then construct a narrow cone to include or exclude charged particles 
reconstructed in the tracking system and/or neutral tt^s reconstructed in the shower 
maximum detector which is inside the electromagnetic calorimeter. 

The question now is: is AR a good choice of cone size definition for hadronic tau 
identification? 

A AR cone has a relativity invariance under a boost along the z-axis. However 
a AR cone does not have geometry invariance. What does a constant AR imply in 
geometry? The top plot of Fig. 15.31 shows three constant isolation annulus at different 
r/ in a uniform r]-(j) space; the bottom plot shows the same three isolation annulus in a 
uniform 6-(f) space after using the function f] = — In tan | to map rj slices to 6 slices. In 
the central region, the isolation annulus is almost unchanged; outside the central region, 
they are severely squeezed, thus Ar/ doesn't have geometry invariance. A(j) doesn't have 
geometry invariance either. Think of one step at the Equator of the Earth and another 
step at the North Pole of the Earth, the former is a tiny one in Ai;^ while the latter is a 
giant one in Acp. A constant AR cone with relativity invariance is not expected to be 
a constant cone with geometry invariance. 

Instead of Et and AR, we can use energy E and three-dimensional separation angle 
a to construct a cone for hadronic tau identification. There are two reasons. 

First, consider a rotation of a solid cone; the geometry invariance of a three- 
dimensional separation angle a is easy to visualize. The unknown boost of the parton- 
parton system along the z-axis doesn't affect the energy measurement of the hadronic 
tau identification at all. Under the known high energy boost, the final particles are 
flying together in a narrow cone. In one case the boost is to the central region, and 
in another case the boost is to somewhere forward or backward. Are these two cones 
geometrically invariant? The answer is yes. 

Second, the correlation of E and a is very strong. The case with the simplest phase 
space of final particles is calculable, see AppendixEl Comparing with a constant cone. 
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Figure 5.3: Lack of geometry invariance in AR cone. 
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a variable cone determined by this correlation can give extra power to suppress the jet 
background for hadronic tau identification. This is described by the "shrinking" cone 
algorithm for hadronic tau identification below. 

5.2.2 The "Shrinking" Cone 

As shown in Fig. l5.2( tau isolation cone, i.e., the outer cone, is a constant 30° (0.525 rad) 
cone. For a particle with definite mass like tau, the bigger the energy, the smaller the 
separation angle of its decay daughters, hence a smaller signal cone which is the inner 
cone in Fig. 15.21 

The tau resonctruction algorithm j25j starts with a seed tower with Et > 6 GeV. It 
adds all of the adjacent shoulder towers with Et > 1 GeV to make a calorimeter cluster. 
The cluster is required to be narrow, i.e., the number of towers < 6. The visible energy, 
denoted as E^s , of the final particles of tau haronic decays is measured by the energy of 
the calorimeter cluster, denoted as -E^/^s^er- Then the algorithm asks a seed track with 
Pt > 4.5 GeV/c to match with the cluster. The matched seed track is a track with the 
highest Pt in the neighbor of the calorimeter cluster. The tau signal cone is constructed 
around the direction of the seed track. The other tracks with > 1 GeV/c, and the 
vr'^s with Pt > I GeV/c which are reconstructed by the strip and wire clusters in the 
CES detector, are included in the tau candidate if they are inside the tau signal cone. 
The size of the tau signal cone is determined by Eyis- 

The phase space of tau hadronic decays is very complicated and the energy depen- 
dence of the signal cone cannot easily be calculated analytically. We use a large MC 
sample of ^ Z —> rr to get this correlation. 

The concept of tau shrinking signal cone at generation level (without underlying 
track or vr*^) is shown in Fig. 15.41 The cone starts out at a constant 10°, and then, if the 
quantity (5 rad) /Ems is less than 10° we use this angle, unless it is less than 50 mrad. 

For reconstructed tracks a cone defined as that shown in Fig. 15.41 is efficient and 
selective against jet backgrounds. However, for n^s, the reconstructed angle can, at 
large visible energies, be larger than 50 mrad. Thus we relax the minimum to 100 mrad. 
With underlying track or tt^, the shrinking cone is shown in the left two plots of Fig. 15.51 
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Figure 5.4: Distribution of maximum angle between tau decay products and tau seed 
track as a function of tau visible decay product energy. The red line indicates the 
half-width of the " shrinking" tau signal cone as a function of energy. 



Inside the tau isolation cone (the outer 0.525 rad cone), the separation angle between the 
farthest track/vr*^ and the seed track is ploted. A tau object between the tau isolation 
cone and the shrinking signal cone is non-isolated and will be removed by isolation cut. 
The right two plots of Fig. 15.5! show how the shrinking cone looks when applied to jets 
reconstructed as tau objects. Comparing with a constant signal cone, the shrinking 
signal cone, a natural consequence of the tau's relativistic boost, dramatically helps to 
reduce jet background in the high mass search. 

5.2.3 Tau Identification Cuts 

Now we can put the seed track in the center of the cone and include in the tau candidate 
all tracks and vr'^s whose direction is within the "shrinking" signal cone. Table l5T] shows 
the list of tau identification cuts using the information about calorimeter cluster, seed 
track, shoulder tracks/vr'^s of the tau candidate. The (tracks + tt^s) threshold is 
not listed because it is not an identification cut and it should be chosen by looking at 
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Figure 5.5: Due to different reconstruction resolutions, the minimum cone sizes of the 
"shrinking" cone for track and vr'' are 0.05 and 0.1 radian, respectively, shown in the 
left two plots for tau. The right two plots show how the "shrinking" cone looks when 
applied to jets reconstructed as tau objects. 



48 



Variable 


Cut 


Note 


Denominator 


\Vdet\ 


<1 


central calorimeter 






9< \zioc\ <230 cm 


fiducial ShowerMax 






>0.2 


electron removal 


s 




>6 GeV/c 


seed track pj- 


10° track isolation 


constant cone 


weaker than shrinking 


DtrklsoWDeg 


m(tracks) 


<1.8 GeV/c2 


weaker than vis. mass 




\zq\ 


<60 cm 


vertex z 




\do\ 


<0.2 cm 


impact prameter 




seed track ax. seg. 


>3x7 


COT axial segments 




seed track st. seg. 


>3x7 


COT stereo segments 




track isolation 


shrinking track cone 


shoulder tracks 




vr*^ isolation 


shrinking tt*^ cone 


shoulder vr'^s 






<2 GeV 


EM cal. isolation 




m(tracks + vr'^s) 


<1.8 GeV/c2 


visible mass 


Numerator 



Table 5.1: Tau identification cuts. 



the trigger cuts applied and by comparing tau identification efficiency with the jet— > r 
misidentification rate. We do not cut on charge because there is an ambiguity in the 
charge for high pT tracks; we do not cut on track multiplicity either because we will 
check track multiplicity to see hadronic tau signature. 

Electron Removal 

Using the requirements discussed above, electrons can be reconstructed as hadronic 
tau objects if they have a narrow calorimeter cluster and a high px seed track. To 
remove electrons we demand that the tau be consistent with having only pions in the 
final state. We define the variable as 

^ = E^I_EeI1) = ^ (5.10) 
p E p 

Fig. 15.61 shows the tau object EM fraction (Eem/E) versus E/p. The top plot is for 
hadronic taus reconstructed as tau objects, and the bottom plot is for electrons recon- 
structed as tau objects. For an ideal hadronic tau and a perfect calorimeter, ^ = 1. For 
an ideal electron, ^ = 0. However, the calorimeter is not perfect and there can be a 
large background from Z ^ ee events. To remove this background we use a very tight 
cut, ^ > 0.2. The remaining background is discussed below. 
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Figure 5.6: Distributions of EM fraction [Eem/E) vs. E/p for hadronic tau and elec- 
tron. ^ > 0.2 is used to remove electron. 

EM Calorimeter Isolation 

The motivation for the EM calorimeter isolation cut is due to vr^ reconstruction 
inefficiency, for example, some CES clusters are not reconstructed as vr'^s if a track 
is nearby. This affects the power of the vr'^ isolation requirement. We add an EM 
calorimeter isolation cut to deal with the remaining jet background. We calculate the 
EM energy in a Ai? = 0.4 cone around the seed track, summing over all EM towers 
which are not members of the tau cluster. Here Ai? is used to calculate the distance 
between the centroid of a calorimeter tower and the seed track because the calorimeter 
tower segmentation is fixed in r/ x space, namely 0.11 x 15° around the central region. 
Since the EM calorimeter isolation cut is strongly correlated with other isolation cuts, 
its marginal distribution is shown in Fig. 15.71 The EM cal. isolation energy versus 
cluster energy plots show that we do not need to use a relative (fractional) cut, which 
is necessary if for high energy tau objects there is significant energy leakage outside tau 
cluster. We instead choose an abosolute cut, £'f™ < 2 GeV. 
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Figure 5.7: Disbutions of EM calorimeter isolation for tau and jet, and distributions of 
EM calorimeter isolation vs. energy of reconstructed tau object for tau and jet. 

Object Uniqueness 

Though not listed in the summary table of tau identification cuts, we note that all 
reconstructed objects in the event are required to be unique. Thus we only apply the 
tau identification cuts to objects not already reconstructed as a photon, electron, or 
muon. In practice, we require that a tau object be 30° away from any identified photon, 
electron, or muon. 

Denominators 

For various subsequent studies presented here we will use specific subsets of the tau 
identification cuts listed in the summary table. The cuts are in cumulative order which is 
important for calculating rates and efficiencies. There are three different denominators 
in Table 15.11 corresponding to three different relative rates, which will be applied on 
different data samples with consistent denominators later. 
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5.2.4 Tau Identification Efficiency 

Table l^^ shows the procedure to measure the tau identification efficiency, using different 
samples. For all of the generated taus, we pick those taus decaying hadronically, and 
consider the central ones in the pseudorapidity range \r]\ < 1 which are able to be 
reconstructed as tau object, called CdfTau in the table. We require the seed track of 
the generated tau to match with the seed track of a reconstructed tau object within 
0.2 radian. Then we apply the tau identification cuts on the reconstructed tau objects 
and calculate tau identification efficiency. 

Fig. 15.81 shows the absolute tau identification efficiency, which includes the effects 
of both reconstruction and identification, vs. tau visible energy, using the Z' sample 
which has a lot of high energy taus. 

5.2.5 Jet— > T Misidentification Rate 

Table ESI shows the procedure to measure the jet— > r misidentification rate, using four 
different jet samples called JET20, JET50, JET70, and JETIOO samples collected with 
different trigger thesholds. The LI tower E^, L2 cluster Et and L3 jet Et trigger 
thresholds in the unit of GeV for a triggered jet in each jet sample are 

• JET20: 5, 15, 20 

• JET50: 5, 40, 50 

• JET70: 10, 60, 70 

• JETIOO: 10, 90, 100 

We use the central jets with |?7| < 1 which may be reconstructed as tau object, called 
CdfTau in the table. We require the central jet to match with a reconstructed tau object 
by requiring that they share the seed tower of the reconstructed tau object. Then we 
apply the tau identification cuts on the reconstructed tau objects and calculate jet— > r 
misidentification rate. 

Fig. 15.91 shows the absolute jet— > r misidentification rate, which includes the effects 
of both reconstruction and identification, vs. jet cluster energy, using JET50 sample. 



52 



Procedure 


W ^TU 


Z TT 


Z' TT 


Denominator 


\^ V L.ilL 


4Q1 ^1 3 

t:?? -LU ±0 


4Q9nnn 






Toil rl O T*I^T1 1 <^ 


OIVOO t 


UO I OOa 










Z ( ooou 




n iL , J 

■'-^absolute 


V CI' LI J.110;uv^J.J. V-^U-X -L CtLl 


86825 




800262 


^CdfTau 


In, 1^1 


oooyy 




7Q77ni^ 




y \zioc\ ^ zou cm 






/ 0O4:UO 






D0o041: 




DDoo^tO 




^ D VorG V / C 


Duyou 


1 1 Q/l c;i 


DOioZo 


10° track isolation 


ououy 


Qfi71 7 

yo ( L( 


04U4oO 


^trklsolODeg 


lllyjjl dbLKh j \ l.O vjrC V / C 




yoooo 
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J^trkMass 
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yoooo 
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\rln\ <r 2 cm 


47975 


93969 


506453 




seed track ax. seg. > 3 x 7 


47822 


93657 


501965 




seed track st. seg. > 3 x 7 


47312 


92666 


494069 




track isolation (shrinking) 


47112 


92042 


475017 




tt'^ isolation (shrinking) 


45687 


89148 


451129 




EfZ < 2 GcV 


43981 


85910 


428641 




m(tracks + 7r°s) < 1.8 GeV/c^ 


43155 


84218 


404105 


Numerator 



Table 5.2: Number of events for tau identification efficiency measurement. 
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Figure 5.8: Tau identification efficiency vs. tau visible energy. 
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Procedure 


JET20 


JET50 


JET70 


JETIOO 


Denominator 


event 


7696880 


1951396 


910618 


1137840 




event goodrun 


4309784 


1213104 


556961 


697231 




iet non-triggered 


21957203 


6071557 


2641643 


2935801 




jet central 


8214991 


2480232 


1127376 


1321840 




jet match CdfTau 


653680 


425086 


189148 


201530 


DcdfTau 


IrirlcA < 1 


643190 


416560 


184996 


196651 




9 < \zioc\ < 230 cm 


611401 


393222 


174474 


184980 




^ > 0.2 


521326 


354504 


159320 


169441 


Dp 


pf?^'^ > 6 GeV/c 


414966 


315384 


145124 


156391 


10° track isolation 


105846 


74425 


36231 


42727 




m(tracks) < 1.8 GeV/c^ 


92475 


63616 


31865 


37709 




\zq\ < 60 cm 


85754 


56951 


28146 


32747 




dol < 0.2 cm 


79889 


51829 


25391 


28994 




seed track ax. seg. > 3 x 7 


78500 


50043 


24293 


27474 




seed track st. seg. > 3 x 7 


71926 


42754 


20058 


21828 




track isolation (shrinking) 


64489 


20679 


7475 


7293 




vr'' isolation (shrinking) 


50886 


13910 


5025 


4965 




< 2 GcV 


41749 


11132 


4073 


3969 




m(( racks + 7r"s) < 1.8 GcV/c- 


3531 1 


79(35 


2879 


2792 


Nuni(;ralor 



Table 5.3: Number of events for jet— > r misidentification rate measurement. 
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Figure 5.9: Jet— > r misidentification rate vs. energy, using JET50 sample 
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Figure 5.10: Discrepancies of jet^ r misidentification rates in JET samples. 



Discrepancies 

To try to minimize trigger bias, we use non-triggered jet only. Based on the LI 
tower Et, L2 cluster Et and L3 jet Et trigger thresholds in each sample, we find all 
of the jets which can satisfy the trigger requirements. The choice of the triggered jets 
in an event in the case of zero, one or more than one jet satisfying trigger requirements 
are 

• If zero, throw away the event 

• If only one, choose that jet 

• If more than one, do not choose any as triggered 

Non-triggered jets are just the jets not chosen as the triggered jet. Even after trying to 
minimize trigger bias by using non-triggered jet only, there are still discrepancies among 
jet— > T misidentification rates obtained from different jet samples, shown in Fig. 15.101 
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Two-Dimensional Parametrization 

There is no doubt that the jet— > r misidentification rate has a very strong depen- 
dence on energy because the tau isolation annulus is a function of energy. To resolve 
the discrepancies among the jet— > r rates, we add another parameter to make a two- 
dimensional parametrization. The second parameter should not be correlated strongly 
with energy, otherwise adding another parameter is meaningless. Given the final par- 
ticles, the transverse size of a jet depends on its boost: jets with a bigger boost have 
smaller size and smaller size jets have higher probability to survive tau identification. 
The relativistic boost 7 is 

7 = - (5.11) 
m 

where E is the energy of the jet which can be measured by its cluster energy in calorime- 
ter, and m is the invariant mass of its final particles. The mass m is not easy to measure 
because some of the final particles can be neutral and leave no track in tracking sys- 
tem. We use cluster mass, which treats each tower in the cluster as a massless photon 
and sums up the photons, as an approximation of m. The cluster mass has a strong 
correlation with energy, while the cluster boost does not. This is shown in Fig. 15.111 
We choose cluster boost as the second parameter. 

In the one-dimensional jet^ r misidentification rate what we see is the average over 
all of the bins of cluster boost. Given the energy of a jet, the average cluster boost is 
different in JET samples, shown in Fig. 15.121 

Now we plot the jet— > r misidentification rate vs. energy, in each boost slice, shown 
in Fig. 15.131 With the new two-dimensional parametrization, the overall discrepancy 
drops down to about 20%. Since the discrepancies are not totally resolved, there are 
other unknown effects. 

5.2.6 Jet— > r Background Estimate 

After applying the full set of tau identification cuts, there will be some jet background 
left because of the huge production rate of jets in pp collisions. The jet— > r misidentifica- 
tion rate and tau identification efficiency are very useful for estimating jet background. 
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Figure 5.11: Distributions of jet cluster: mass vs. energy, and boost vs. energy, in 
JET50 sample. 
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Figure 5.12: Profiles of jet cluster boost vs. cluster energy in JET samples. 
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Figure 5.13: Jet— > r misidentification rate vs. energy in JET samples and in jet cluster 
boost slices. 
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To estimate the jet background, the starting point is not jets, or tau candiates, but 
tau candidates with at least electron removal, with a very tight ^ > 0.2 cut apphed. 
Muons usually cannot have enough energy to make a tau cluster in the calorimeter. We 
have two general equations, 

Before full tau ID: N = + iV^'^* (5.12) 
After full tau ID: N = N'' + iV^'^* = eAT^ + /iV^''* (5.13) 

where / is jet— > r misidentification rate and e is tau identification efficiency. Both are 
relative in a sense that they are relative to the starting point chosen as "Before full tau 
ID". The solution is 

N^^^ = — ^(eiV - N). (5.14) 

Fig. 15.1^ is a demonstration of picking one bin and using the formula to estimate jet 
background. This is only an example because the parametrization of the relative rates 
is a one-dimensional function of energy. For the jet— > r misidentification rate there is 
a better parametrization, i.e., the two-dimensional function of energy and boost. 

Implementation 

The actual implementation is done on an event-by-event basis. For a tau object in 
an event under consideration, the knowns are: N = 1, e, f and whether this tau object 
passes the full set of the tau identification cuts. If it does, = 1; otherwise, = 0. 
For the two cases, the weight to be a jet is estimated as 

f 

If not passing the fuh tau ID cuts: w^^* = — (e - 0) (5.15) 

e - / 
f 

If passing the fuh tau ID cuts: u^"^ = ^!—-{e - 1) (5.16) 

^ J 

In terms of coding, it means the rest of full tau identification cuts are replaced by the 
weight w-^'^*. We sum up the weights of all the events in the sample, and get the jet 
background estimate N^^^, 
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Special Case 

This method actually needs both the jet— > r misidentification rate / and the tau 
identification efficiency e. The main idea is to remove the contribution from any real 
tau signal in jet background estimate. 

The special case is that if we start with a jet-dominated sample and / is much 
smaller than e, then we can suppress signal by replacing tau identification cuts with 
the jet— > r misidentification rate, 

iVi^* = /iVJe* ^ fN (/ < e) (5.18) 
5.3 Tau Scale Factor Using W ^ tu 

In this section, we apply tau identification cuts to select hadronic taus va. W ^ tv 
events, estimate jet— > r misidentification background, study tau identification scale 
factor and compare tau distributions in data and MC simulation. 

5.3.1 Data/MC Scale Factor 

The scale factor for a set of cuts quantifies and corrects for the difference between 
data and MC simulation. It should be multiplied on MC to get the scaled efficiency 
consistent with the efficiency in data. Fig. 15. I.'tI shows lepton flow in data and in MC, 
and lepton data/MC scale factors. 

Ratio of Efficiencies 

A data/MC scale factor is defined as the ratio of efficiencies, 

fdata/MC = (5.19) 

where eMC is the efficiency in MC which is straightforward to obtain because the MC 
simulation has the true information of particle identity, and edata is the efficiency in 
data, which can be a challenge to measure. 

In the electron or muon case, we can use electron or muon pairs from the Z boson 
peak, which gives us a pure sample with negligible background in real data. This 
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Figure 5.15: Lepton in data and in MC, and lepton data/MC scale factors. 



is so reliable that we can use it as "standard candle" to calibrate detector and even 
measure luminosity. We select one leg to satisfy the trigger requirements in data, and 
ask whether the second leg passes the set of cuts, and thereby get the efficiency in data. 

Ratio of Numbers 

Due to the missing energy from the neutrino in tau decays, the tau pair mass at the 
Z boson peak is severely broadened. Instead, we will use W ^ tv to select a relatively 
clean tau sample. There is no second leg to get efficiency data/MC. We use the method 
of absolute number data/MC, 

f data/MC = (5.20) 

where umc is the absolute number olW ^ tv events in MC normalized to the lumi- 
nosity of data, and Udata is the number of W ^ tv events observed in the data after 
subtracting backgrounds. 

5.3.2 W -^Tu Selection 

We select W ^ tv events by using a data sample from the TAU_MET trigger which 
requires: 
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Figure 5.16: Distributions of offline luminosity vs. good run sequence and L3 cross 
section vs. good run sequence, in the data sample from TAU_MET trigger. 

• level 1 trigger (LI) > 25 GeV 

• level 3 trigger (L3) tau Et > 20 GeV 

where (a) LI is based on a tower threshold of 1 GeV for a fast calculation; (b) for 
L3 tau, the cuts |%ei| < 1) 10° track isolation and m(tracks) < 2 GeV/c^ are applied 
in the trigger. 

The top plot of Fig. 15.161 shows that the integrated luminosity of the good runs is 
72 lb 4 pb~^, and the bottom plot shows the L3 cross section is reasonablly flat (no 
sudden drop to zero), thus all of the good runs are present in the data file. 

The offline selection cuts are: 

• Mono jet 

• > 30 GeV 

• Tau pT(tracks + 7r°s) > 25 GeV/c 

where (a) monojet selection requires one central cone 0.7 jet with Ir/^etl < 1 and Et > 
25 GeV, no other jets with Et > 5 GeV anywhere; (b) offline ^t is obtained from the 



Correc 


ti< 
7 

i 


an 

2 1 

i 


+i.< 

4 p 


}%, ern 


or 6% 


Enl 
Intc 


ries 

;gral 7( 


460 
).73 


TAU_MET 

■ 1 1 ■ 1 ■ 1 L. 


■ 1 ■ 


/\ 

—1 ■ ■ i_ 


V|A 


. il.lll 1 
—1 ■ 1 1 


ll, 



63 



vector sum of Et for towers with Ej- > 0.1 GeV; (c) in addition to tau pT threshold, 
the whole set of tau identification cuts under study will be applied on the offline tau 
candidates. 

The monojet cut dramatically helps clean up the data sample. But, to get the 
estimated umc of W ^ tu events, we need to study the monojet cut and the LI > 
25 GeV trigger efficiency for monojet-type events. 

Monojet Selection 

The monojet selection essentially requires there is no other underlying jet with Et > 
5 GeV. We select Z ^ fifi events, count the number of cone 0.7 jets with E^ > 5 GeV, 
no r] cut, and 0.7 radian in AR away from muons. 

The Z /i/x selection cuts are: (a) cosmic veto ^], (b) one tight muon and one 
track with px > 20 GeV/c, (c) opposite charges, (d) track \zq{1) — zq{2)\ < 4 cm, and 
(e) 80 < < 100 GeV/c^. We require one tight muon and one track, instead of two 
tight muons to get higher statistics. The track is required to be of minimum ionisation 
particle (MIP) type. Both the tight muon and the track requires tau-like track isolation 
which is to mimic the isolated tau va. W ^ tv events. 

We use a data sample from a trigger designed to select "muon plus track" events 
which have // with pT > 8 GeV/c plus another charged track with pT > 5 GeV/c. 
We select 5799 events with negligible background which is confirmed by the negligible 
number of same-charge muon pair events. There are 2152 events in the zero jet bin. 
The fraction of zero jet events in the data is 2152/5799 = 0.371. 

We use about 500K MC events. 46297 events survived after the same selection cuts 
as in data. There are 20149 events in the zero jet bin. The fraction of zero jet events 
in the MC is 20149/46297 = 0.435. 

The number of jets distribution in data and in MC are shown in Fig. 15.171 So 
W ^ TV monojet data/MC scale factor is 

_ 2152/5799 _ 0.371 2l\ 
hata/MC - 20149/46297 " 0435 " ^'^^ ^ 

The uncertainty is statistical only. 
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LI $T > 25 GeV 

The TAU_MET trigger triggers directly on tau objects, and so there is no marginal 
trigger efficiency from TAU side. But there is marginal trigger efficiency from MET side: 
LI uses a 1 GeV tower threshold, and offline uses a 0.1 GeV tower threshold. 

We use JET20 data to study this trigger efficiency. The event topology is monojet- 
like, since here that is what we are interested in. The LI ^j- > 25 GeV trigger efficiency 
vs offline for monojet event is shown in Fig. 15.181 It is a slow turn-on due to a large 
tower threshold. An offline > 30 GeV cut is not fully efflcient. 

5.3.3 Tau Scale Factor 

After all of the above, we count the absolute number W ^ riy events n^iata aud umc 
for total integrated luminosity 72 pb~^. Their ratio will be the tau scale factor. 

• To get Tiiiata-, we will use the data sample from the TAU_MET trigger. We apply 
the offline cuts to get the observed number oiW^rv candidates, and subtract 
various backgrounds. 

• To get Umc-, we will use W ^ tu MC simulation. We apply the offline cuts, 
multiply the number of accepted events by the monojet scale factor and the 
trigger efficiency, and normalize to 72 pb~^. 

The main source of backgrounds are W — > eu, W fiv, Z/^* — > rr, and jet 
background. 

We will use MC simulation to get W — > ez^, W — > f^u, and Z/^* tt backgrounds. 
We apply the offline cuts, multiply the number of accepted events by the monojet scale 
factor and the trigger efflciency, and normalize to 72 pb~^. For the normalization in MC, 
a-B{W lu) is 2700 pb 20 , and a ■ B{Z /-f* II) is 326 pb with m^/^. > 30 GeV/c^, 
which is obtained from the measured value 250 pb j20j at the Z boson mass peak with 
66 < mz/^* < 116 GeV/c^ and normalizing the Z/^* —>■ II mass spectrum generated 
by PYTHIA. 
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Figure 5.17: Distributions of the number of jets in Z ^ fi/j, data and MC. 
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Figure 5.18: LI > 25 GeV trigger efficiency vs. offline for monojet event. 
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W ^TU 


W ^eu 


W ^ fxu 


Z/Y TT 


etau08 




signal 


bkgd 


bkgd 


bkgd 


jet bkgd 


observed 


event 


491513 


1480550 


760457 


492000 


3747680 


hadronic tau 


319517 


N/A 


N/A 


N/A 


TAU_MET 342164 


mono jet 


11368 


192806 


7557 


7311 


23818 


$T > 30 GeV 


4874 


154256 


4535 


3149 


17490 


tau ID 


1982 


319 


130 


1230 


DtrkMass 1519 


monojet SF 


1684.7 


271.2 


110.5 


1045.5 




tau ID 


trigger eff. 


1622.1 


267.0 


107.6 


1012.3 


81.8 


814 


normalized 


638.8 


34.9 


27.4 


48.3 


81.8 


814 



Table 5.4: Expected number of events for the signal, backgrounds and observed number 
of W ^ TV events. 



The jet background will be estimated directly from the data by applying the relative 
jet— > r misidentification rate and the relative tau identification efficiency. Since the 
cuts \r](iet\ < I7 10° track isolation and m(tracks) < 2 GeV/c^ are applied in the trigger, 
we use the relative rates up to the denominator DtrkMass- Then we just follow the 
implementation described in section [5.2.61 

Table 1231 shows the procedure to estimate the contributions from signal and back- 
grounds estimated from MC, the jet background estimated from data, and the observed 
number of events in data. 

The uncertainties include 

• statistical uncertainty, 

• monojet scale factor: 2%, 

• luminosity: 6% ^27j, and 

• a ■ B{W — > Iv) and a ■ B{Z/^* II), 2%, aside from luminosity uncertainty |20j . 

Since there are discrepancies among the jet— > r misidentification rates obtained 
from different jet samples, we use the average jet^ r misidentification rate to get a 
central value of 81.8 events. The estimates using the individual jet^ r misidentification 
rate from JET20, JET50, JET70, and JETIOO samples are 90.3, 67.1, 72.8, and 66.3, 
respectively. We take the biggest difference as the the uncertainty for jet background: 
|(66.3-81.8)/81.8| = 18.9%. 
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W ^TU 

W ^ev 
W ^ fiu 

Z/-f* TT 

Jet— > T 


638.8 ±42.7 
34.9 ±2.9 
27.4 ± 3.0 
48.3 ±3.3 

81.8 ± 15.5 


Expected 


831.2 ±45.7 


Observed 


814 



Table 5.5: Cross check for the numbers of W ^ ti' events. 



The numbers and the uncertainties of each channel are summarized in Table 15.51 
We now arrive at the tau scale factor as follows: 



data/MC 



riMC 

nobs. — nwZ bgs — njet bg 

n 



sig. 

814 - (34.9 ± 27.4 ± 48.3) - 81.8 
638.8 

= 0.97 ±0.10 (5.22) 

with statistical uncertainty and all of the systematic uncertainties. 

Lastly we put signal and background together, and show the W ^ tv kinematic 
distributions in data and MC in Fig. 15.191 The agreement between data and MC is 
very good. 

5.4 Electron Identification 

Identification of electrons is based on the energy it deposits in the calorimeter, its 
track in the COT, and its position in the CES. The central electron reconstruction 
algorithm [2Hj starts with clusters in the CEM detector. The electromagnetic towers 
are ordered in Et and the highest Et tower that has not yet been clustered is taken 
as a seed. The available shoulder towers are added to the cluster if they are adjacent 
in r/ to the seed, and the clusters are restricted to two towers. The default threshods 
for seed towers and shoulder towers are 3.0 and 0.1 GeV, respectively. For the leading 
electrons used in our analysis with e + t^ channel, the threshods for seed towers and 
shoulder towers are set to be 8 and 7.5 GeV, respectively. Then we associate tracks 



68 



400 



(0 

c 

0) 

> 

UJ 



200 



CDF Run II 72 pb"* 
Preliminary 
Tau Identification 



QW^TV 
□W->ev 

□ Z^TT 

■Jet^T 



> 
O 

CM 
(0 

c 

0) 

> 

UJ 




100 



30 40 50 60 

p-rCtracks + -nPs) [GeV/c] 





^200 



Efs-S [GeV] 



-1 1 
Charge 



Figure 5.19: Distributions of hadronic tau identification using W ^ tv events for data 
(points) and predicted backgrounds (histograms). 
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Variable 


Cut 


Note 


Denominator 


region 


==0 


CEM 




fiducial 


==1 


fiducial Xces, Zqes 




ko 


<60 cm 


vertex z 


Probe 


track ax. seg. 


>3x7 


COT axial segments 




track st. seg. 


>3x7 


COT stereo segments 




cal. isolation 


<0.1 


cone 0.4 




Ehad/ Eem 


<0.055+0.00045x£; 


had. /em. 




E/p 


<4 (for Et <100 GeV) 


cal. /track with brem. 




^shr 


<0.2 


lateral shower profile 




\AX\ 


<3 cm 


Xtrack — Xces 




\AZ\ 


<5 cm 


Ztrack — ZcES 




conversion veto 


|AXy|<0.2 cm, and 


separation, and 






Acot6l|<0.04 


parallel 


Numerator 



Table 5.6: Electron identification cuts. 



with the candidate cluster. For all of the tracks associated, the one with highest px is 
chosen as the matched one. The CES strip and wire clusters are associated with the 
CEM cluster if they are reconstructed in the same wedge. The "best-matching" CES 
cluster is the one seeded by the matched track. 



5.4.1 Electron Identification Cuts 

The electron identification j^H] cuts, and the conversion veto i.30| cuts to remove elec- 
trons from photon conversion, are listed in Table 15.61 The Et and px thresholds are 
not listed because they depend on the process and trigger sample. The probe electron 
must be a fiducial CEM electron and pass the vertex z cut. 



5.4.2 Electron Scale Factor 

The electron identification scale factor is the ratio of the efficiency in data/MC. The 
data sample is from the TAU_ELE trigger which requires an electron with Et > 8 GeV, 
> 8 GeV/c and an isolated track with pT > 5 GeV/c. We study the electron scale 
factor versus Et- 

• For medium-ET (between 5 GeV and 20 GeV) electrons, the MC uses electrons 
from Z — > rr — > eX, and in the real data we use the second leg after selecting 
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T ee. We require the probe electrons have Et > 5 GeV and pr > 5 GeV/c in 
both the real data and the MC. 

• For high-£^r (above 20 GeV) electron, the MC uses electrons from Z ee, and 
in the real data we use the second leg after selecting Z ee. We require the 
probe electrons have Et > 20 GeV and pT > 10 GeV/c in both the real data and 
the MC. 

The procedure to select T — > ee events is: 

• Require a tight electron with Et > 8 GeV, pt > 8 GeV/c which are the trigger 
requirements and the electron identification cuts. 

• Require a probe electron with Et > 5 GeV, pT > 5 GeV/c. 

• Same-sign pair will be used later for fitting the slope of background and opposite- 
sign pair will be used later for fitting signal -|- background. 

• Require the invariance mass of the ee pair to lie in the range (0, 20) GeV/c^. 
The procedure to select Z ^ ee events is: 

• Require a tight electron with Et > 20 GeV, pT > 10 GeV/c and the electron 
identification cuts. 

• Require a probe electron with Et > 20 GeV, pT > 10 GeV/c. 

• Require opposite sign. 

• Require the invariance mass of the ee pair to lie in the range (75, 105) GeV/c^. 
The procedure to select the second leg is: 

• Require exactly one T or Z boson. 

• If there is one tight electron, the probe electron is the second leg. 

• If there are two tight electrons, both are used as second leg. 
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data 




MC 






T — > ee 


Z — > ee 




Z — > T^Tx 


Z — > ee 




event 


11922805 


event 


492000 


398665 




good run 


9103020 










triggered 


5575584 










unique 


5310963 










process 


10373 


4534 


electron 


175515 


797330 




second leg 


10770 


7973 


match 


36733 


204680 


Probe 


track ax. seg. > 3 x 7 


10687 


7946 


same 


36665 


204266 




track st. seg. > 3 x 7 


10165 


7721 


same 


36448 


202907 




cal. isolation < 0.1 


2797 


7484 


same 


32094 


197065 




Ehad/Eem < 0.055+0.00045S 


2553 


7427 


same 


31290 


194746 




E/p < 4 (for Et < 100 GeV) 


2551 


7379 


same 


31187 


194133 




Lshr < 0.2 


2331 


7318 


same 


30240 


188653 




\AX\ < 3 cm 


2304 


7198 


same 


30028 


186360 




AZ < 5 cm 


2292 


7189 


same 


29976 


186222 




conversion veto 


2249 


6878 


same 


29714 


181449 


Id 



Table 5.7: Number of events for electron identification efficiency measurement. 



Then we apply the set of electron identification cuts under study on the second leg 
electrons in data, and on the probe electrons in the MC. The result of the procedure is 
shown in Table 1^771 

For the Z — > ee selection in the real data, the backgrounds in the sample with a 
tight electron plus a probe electron and in the sample with two tight electrons are both 
negligible which is confirmed by the negligible number of same-sign events in these two 
samples. 

For the T — > ee in the real data, the backgrounds in the sample with tight electron 
plus probe electron and in the sample with two tight electrons are both significant. 
The same-sign samples provide the shapes of the invariant mass distribution of the 
backgrounds, which are taken as the slopes of linear backgrounds. Then in the opposite 
sign samples we fit the invariant mass distributions by the "Crystal Ball" function 
plus a linear background. The T — > ee invariant mass distribution has a Bremsstrahlung 
tail at lower mass side where at least one of the electrons radiates. The "Crystal Ball" 
line-shape serves to model this Gaussian core with a power-law tail. The yield of signal 
is obtained by the entries in the histogram subtracted by the integral of the linear 
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background. 

Up to this point, all of the T ^ ee candidates in the mass window (0, 20) GeV/c^ 
are accepted. We then subtract background, as shown in Fig. I5.2UI The plot only 
shows the mass window (4, 15.5) GeV/c^. The fit is performed in the mass window 
(5, 12.2) GeV/c^. The fitting result is N(e + probe) = 818.0, N(e + Id) = 644.4, 
efficiency = 78.8%. 

Now we put everything together to get the electron scale factor vs. Et and perform 
a fit in Et- This is shown in Fig. 15.211 Data: the medium Et electrons (5, 20) GeV 
are from the second leg of T — > ee; the high Et electrons (30, 100) GeV are from the 
second leg of Z ^ ee; there is a gap (20, 30) GeV which has very low statistics and is 
not used. MC: the medium Et electrons (5, 20) GeV are from the probe electrons of 
Z ^ TT ^ eX; the high Et electrons (30, 100) GeV are from the probe electrons of 
Z ee. In each Et bin, the efficiency in data divided by the efficiency in MC gives 
scale factor in that Et bin. For all of the Et bins, the scale factor is flat. A fit by a 
polynomial of degree 0, which is exactly the same as the weighted average, gives a scale 
factor 0.974 ± 0.004. 

There are two bins (45, 50) GeV and (50, 100) GeV with efficiency close to 100%, 
in data and MC. The binomial uncertainty in this case is always close to zero and 
underestimated. This propagates to the scale factors in those two Et bins, and finally 
propagates to the weighted average. There is also uncertainty in the (5, 20) GeV bin 
due to T ^ ee background subtraction. This uncertainty is not estimated. We assign 
a conservative 4% uncertainty for electron scale factor j32j : 

/1WMC = 0-97 ±0.04 (5.23) 

5.5 Muon Identification 

Muon reconstruction uses information from tracking, the calorimeters and the muon 
chambers. The momentum is measured by the curvature of the muon trajectory bent 
by magnetic field in tracking system. Muons behave as minimum ionizing particles and 
they are the only charged particle that can travel through the large amount of material 
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m [GeV/c1 



m [GeV/cl 



Figure 5.20: Distributions of the invariant mass of T — > ee for medium Et electron 
identification efficiency measurement in data. Same-sign samples provide the slopes 
of the linear backgrounds. The numbers of T ^ ee signal events are obtained from 
fitting the histograms with the "Crystal Ball" function plus a linear background in the 
oppisite-sign samples. 
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Figure 5.21: Electron scale factor vs. Et- This is obtained from dividing the efficiency 
in data by the efficiency in MC. 



74 



in calorimeter with a very small energy loss. Muons are not stable, but they are so long 
lived that they can reach the muon chamber, leave hits there, and continue to travel 
and decay outside the detector. These features allow a rather simple and clean muon 
identification. 

5.5.1 Muon Identification Cuts 

The muon identification cuts are listed in Table IKl^ We use COT-only tracks and 
add the beam constraint to the track. The Et and pT thresholds are not listed because 
they depend on the process and trigger. For data/MC scale factor studies, we require 
the track to be fiducial which means that the track is headed in a direction that will 
lead it to hit enough chambers for a stub to be reconstructed. All three subdetectors 
CMU, CMP, and CMX which are used in this analysis require 3 hits in 3 different layers 
for a stub to be reconstructed. And we will study two kinds of data/MC scale factors: 

• Muon identification scale factor. A fiducial stub muon and vertex z cut are 
required for the probe muon for this study, called "Probe (Id)" in Table EHl 

• Marginal muon reconstruction scale factor. We require a fiducial track and a 
stubless muon, which also has the information of energy loss in calorimeter, for the 
probe muon for this study, called "Probe (Rec)" in Table It is not necessary 
to have hits in the muon chambers. The vertex z cut, calorimeter isolation cut, 
EM energy cut, and hadronic energy cut are required. Then we check if this track 
has a muon stub. The default track px threshold to make a stubless muon is 10 
GeV; we lower it to 5 GeV to allow more medium pT stubless muons. 

5.5.2 Muon Scale Factor 

The muon identification scale factor is the ratio of the identification efficiency in real 
data to that in MC. The muon marginal reconstruction scale factor is the ratio of the 
marginal reconstruction efficiency in data/MC. The data sample is from the TAU_CMU 
trigger which requires a CMUP muon with pT > 8 GeV/c and an isolated track with 
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Variable 


Cut 


Note 


Probe 




<60 cm 




vertex z 


Probe 


(Id) 


cal. isolation 


<0.1 




cone 0.4 








<2+max(0,(p-100)x 


0.0115) GeV 


EM energy 






Ehad 


<6+max(0,(p-100)x 


0.028) GeV 


had. energy 


Probe 


free.) 




<0.2 cm 




impact parameter 






track ax. seg. 


>3x7 




COT ax. seg. 






track st. seg. 


>3x7 




COT St. seg. 






I^^CMUI 


<3 cm (for CMUP) 




^track - ^CMU 






lAxcMpI 


<5 cm (for CMUP) 




^track - ^CMP 








<6 cm (for CMX) 




^track - ^CMX 






PCOT 


>140 cm (for CMX) 




COT exit radius 







Table 5.8: Muon identification cuts. 



PT > 5 GeV/c. (A CMUP muon is required to have stubs in both CMU and CMP). 
We study the muon scale factors versus muon pT- 

• For medium pT (between 5 and 20 GeV/c) muons, the MC uses muon from 
Z ^ TT ^ fj,X, data uses the second leg after selecting T fifi. We require the 
probe muons have px > 5 GeV/c in both the real data and the MC. 

• For high px (above 20 GeV/c) muons, the MC uses muon from Z /i//, and for 
the data we use the second leg after selecting Z fifi. We require the probe 
muons have px > 20 GeV/c in both the real data and the MC. 

The procedure to select T — > /i/z events is: 

• Cosmic veto ^H] . 

• Require a tight CMUP muon with px > 8 GeV/c which are trigger requirements 
and the CMUP muon identification cuts. 

• Require a probe muon with pT > 5 GeV/c. 

• Require \zo{l) — zq{2)\ < 4 cm. 

• Require opposite sign. 

• Mass window (7, 13) GeV/c^. We will use side band for background subtraction. 
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The procedure to select Z ^ ^jj. events is: 

• Cosmic veto. 

• Require a tight CMUP muon with px > 20 GeV/c and the CMUP muon identi- 
fication cuts. 

• Require a probe muon with pT > 20 GeV/c. 

• Require |2;o(l) — zq{2)\ < 4 cm. 

• Require opposite sign. The neghgible number of same sign events confirms that 
background is neghgible. 

• Mass window (80, 100) GeV/c^. 
The procedure to select the second leg is: 

• Require exactly one T or Z boson. 

• If one tight muon, the probe muon is the second leg. 

• If two tight muons, both are used as second leg. 

Muon Identification Scale Factor 

In the muon identification scale factor study, we apply the set of muon identification 
cuts under study on the second leg muons in data, and on the probe muons in the MC. 
Table shows the procedure in data and Table I^TTUl shows the procedure in MC. 

Up to this point, all of the T /u/x candidates in mass window (7, 13) GeV/c^ are 
accepted. Now we break the probe into two px bins 5 < < 10 GeV/c and 10 < px < 
20 GeV/c. Fig. 15.221 shows the distributions of the pair mass of the first leg and the 
second leg in each px bin of the second leg, for CMUP probe. We see three clear peaks 
at about 9.5, 10 and 10.3 GeV/c^. This is the signature of T ^ fi/j,. Now we subtract 
the linear background. The signal mass window is defined as (9.2, 10.6) GeV/c^. We 
use a side-band method: yield = entries in (9.2, 10.6) — entries in (7.8, 8.5) — entries in 
(11.3, 12) GeV/c^ mass windows. In 5 < pT < W (10 < pr < 20) GeV/c bin, N(muon 
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nil 


Z ^ 


nil 




event 


11922805 




good run 




9103020 






triggered 




1881529 






unique 




1800059 






process 


971 


2025 




second leg 


1047 


2805 






CMUP 


CMX 


CMUP 


CMX 






762 


285 


1820 


985 


Probe 


\do\ < 0.2 cm 


758 


283 


1816 


982 




track ax. seg. > 3x7 


756 


283 


1814 


978 




track st. scg. > 3x7 


739 


280 


1777 


955 




cal. isolation < 0.1 


527 


194 


1750 


947 




Eem < 2+max(0,(p-100)x0.0115) GeV 


525 


191 


1700 


929 




Ehad < 6+max(0, (p- 100) X 0.028) GeV 


524 


191 


1670 


907 




I^^CMUI < 3 cm (lAxQjvixl < 5 cm) 


427 


129 


1590 


877 




^^CMP ^ ^ ^™ (/'COT ^ *^™) 


287 


113 


1560 


753 


Id 



Table 5.9: Number of events for muon identification efficiency measurement in data. 





Z 




Z 


nil 




event 


492000 


405291 




muon 


170596 


810582 




match 


32388 


170382 






CMUP 


CMX 


CMUP 


CMX 






20516 


11872 


107481 


62901 


Probe 


\do\ < 0.2 cm 


20499 


11827 


107410 


62692 




track ax. seg. > 3x7 


20490 


11732 


107348 


62184 




track st. seg. > 3x7 


20417 


11601 


106985 


61463 




cal. isolation < 0.1 


18361 


10448 


104247 


59929 




Eem < 2+max(0,(j?-100)x0.0115) GeV 


18055 


10261 


100115 


57680 




Ehad < 6+max(0, (p- 100) X 0.028) GeV 


17871 


10080 


97906 


55799 




I^^CMUI < 3 cm (|A,TQ]Y/Ixl < 5 cm) 


16799 


8994 


97707 


55617 




I^^^CMPl < 6 cm (pcOT > 1^0 cm) 


14198 


7419 


96788 


46059 


Id 



Table 5.10: Number of events for muon identification efficiency measurement in MC. 
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T ^ 




Z 


nil 




event 


11922805 


goodrun 




9103020 






triggered 




1881529 






unique 




1800059 






process 


691 


1861 


Probe 


second leg 


760 


2583 


CMUP 


CMX 


CMUP 


CMX 


570 


190 


1718 


865 


has stub 


474 


170 


1569 


838 


Rec. 



Table 5.11: Number of events for muon reconstruction efficiency measurement in data. 

+ probe) = 410 (85), N(muon + Id) = 168 (56), we get efficiency = 41.0% (65.9%). 
We put everything together to get the CMUP muon identification scale factor vs. pT 
and perfom a fit in p^- This is shown in Fig. 15.231 

Fig. 15.241 shows the mass distribution of muon pair in each px bin of the second leg, 
for CMX probe. In 5 < < 10 (10 < pr < 20) GeV/c bin, N(muon + probe) = 126 
(32), N(muon + Id) = 57 (19), we get efficiency = 45.2% (59.4%). Fig. E^S] shows the 
procedure to get the CMX muon identification scale factor. 

Analogous to the electron scale factor study, we assign a conservative uncertainty 
of 4%. The resulting identification scale factors are 0.93 it 0.04 for CMUP muon, and 
1.03 ± 0.04 for CMX muon. 

Muon Reconstruction Scale Factor 

In the reconstruction scale factor study, the probe is a stubless muon which may 
or may not have a stub in the muon chamber associated with the fiducial track. It 
must have passed the vertex z cut, calorimeter isolation, EM energy cut, and hadronic 
energy cut. For such a probe, we check if it has a stub. Table IF!TT] shows the procedure 
in data and Table I^TT^ shows the procedure in MC. 

We break the second leg muon into two pT bins: 5 < p-r < 10 GeV/c and 10 < < 
20 GeV/c. Fig. 15.261 shows the distributions of the pair mass of the first leg and the 
second leg in each pT bin of the second leg, for CMUP probe. We use the side-band 
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Figure 5.22: Distributions of the invariant mass of T ^ /i/x for medium pj- CMUP 
muon identification efficiency measurement in data. The three peaks are signature of 
T fifx. The left two plots are for CMUP muons with 5 < pT < W GeV/c. The right 
two plots are for CMUP muons with 10 < px < 20 GeV/c. Side-band method is used 
for background subtractions. 
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Figure 5.23: CMUP muon identification scale factor vs. pT- This is obtained from 
dividing the efficiency in data by the efficiency in MC. 
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Figure 5.24: Distributions of the invariant mass of T ^ /u/x for medium px CMX muon 
identification efficiency measurement in data. The three peaks are signature of T ^ ^u/x. 
The left two plots are for CMX muons with 5 < pj- < 10 GeV/c. The right two plots are 
for CMX muons with W < px < 20 GeV/c. Side-band method is used for background 
subtractions. 
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Figure 5.25: CMX muon identification scale factor vs. pT- This is obtained from 
dividing the efficiency in data by the efficiency in MC. 
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Z ^ fin 


Probe 


event 


492000 


405291 


muon 


170596 


810582 


match 


27109 


149411 


CMUP 


CMX 


CMUP 


CMX 


17334 


9775 


95503 


53908 


has stub 


16672 


9709 


93044 


53827 


Rec. 



Table 5.12: Number of events for muon reconstruction efficiency measurement in MC. 

method to do background subtraction. In 5 < < 10 (10 < pT < 20) GeV/c bin, 
N(muon + probe) = 307 (65), N(muon + stub) = 272 (59), we get efficiency = 88.6% 
(90.8%). We put everything together to get the muon reconstruction scale factor vs. 
PT and perform a fit in pT- This is shown in Fig. 15.271 

Fig. 15.281 shows the mass distribution of muon pair in each pT bin of the second leg, 
for CMX probe. In 5 < pT < 10 (10 < pr < 20) GeV/c bin, N(muon + probe) = 92 
(22), N(muon + stub) = 85 (21), we get efficiency = 92.4% (95.5%). Fig. ESi shows 
the procedure to get the CMX muon reconstruction scale factor. 

As in the electron scale factor study, we assign a conservative systematic uncertainty 
of 4%. The results of the reconstruction scale factors are 0.94 it 0.04 for CMUP muon, 
and 0.97 ± 0.04 for CMX muon. 

We summarize the muon reconstruction and identification scale factors with uncer- 
tainties: 



fCMUP rec 
Jdata/MC 

fCMUP Id 
Jdata/MC 

fCMX rec 
Jdata/MC 

fCMX Id 
Jdata/MC 



0.94 ± 0.04 

0.93 ±0.04 
0.97 ±0.04 
1.03 ±0.04 



(5.24) 
(5.25) 
(5.26) 
(5.27) 



5.6 Missing Transverse Energy 

Weakly interacting particles such as neutrinos of the Standard Model and the light- 
est super symmetric particle (LSP) predicted in new physics, deposit no energy in the 
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Figure 5.26: Distributions of the invariant mass of T ^ /i/x for medium pT CMUP 
muon reconstruction efficiency measurement in data. The three peaks are signature of 
T fifi. The left two plots are for CMUP muons with 5 < pT < W GeV/c. The right 
two plots are for CMUP muons with 10 < px < 20 GeV/c. Side-band method is used 
for background subtractions. 




20000 



20 40 60 80 100 
Pt [GeV/c] 



10000 
1 

0.5 
0, 



■ Z->TT 






■n. 


■ 1 ■ ■ ■ 1 ■ ■ ■ 1 ■ ■ ■ 

/] 

r 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 


- H,4-»-l ►•+«+»+»H • 

... 1 ... 1 ... 1 ... 1 .. . 



20 40 60 80 100 
Pt [GeV/c] 



0) 0.5 

(0 

o 
V) 









l-J 


-j^«--*-«-|^-»-r 






pO 0.939+0.00641 




CMUP Rec SF = ^.''^'^ 


1 


"=MC 

1 ... 1 ... 1 .. . 



20 



40 60 
Pt [GeV/c] 



80 



100 



Figure 5.27: CMUP muon reconstruction scale factor vs. pT- This is obtained from 
dividing the efficiency in data by the efficiency in MC. 
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Figure 5.28: Distributions of the invariant mass of T ^ fi/j, for medium px CMX 
muon reconstruction efficiency measurement in data. The three peaks are signature of 
T fj,fi. The left two plots are for CMX muons with 5 < pT < W GeV/c. The right 
two plots are for CMX muons with 10 < pT < 20 GeV/c. Side-band method is used for 
background subtractions. 
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Figure 5.29: CMX muon reconstruction scale factor vs. pT- This is obtained from 
dividing the efficiency in data by the efficiency in MC. 
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calorimeters. Minimum ionizing particles such as muons leave little energy in the 
calorimeters. When present these cause a significant vector sum of the transverse 
energy of all of the detected particles. The imbalance, i.e., the negative of the vector 
sum in the transverse plane corresponds to the missing transverse energy (^t)- 

Since I^t measures the vector sum of all of the momentum of particles escaping 
detection in the calorimeters, there is no information on the energy and direction of an 
individual particle or how many particles escaped detection. With many such particles 
in an event there is also a chance that their transverse momenta cancel each other. 

There is an instrumental source of because the calorimeters are not perfect. 
There are crack regions due to the support structures, and the transition regions be- 
tween components, for example from the central calorimeters to the plug calorimeters. 
The probability that all the energy of a particle is undetected is rather small. But 
QCD processes have a large production rate. Some of the jets can have a lot of energy 
undetected and make a significant ^t- 

In our high- mass tau pair analysis, we will use an IpT cut and several other kinematic 
cuts related to IpT- To get the uncertainty due to the instrumental IpT, we should get 
the distributions in data and MC, and compare the same variable. 

In the real data, the physics processes Z — > ee and 7 -|- jet, which have zero true 
missing energy, can be used to study the effect of the instrumental IpT- The latter is a 
better choice for our purpose because hadronic taus in the calorimeters are more like 
jets than electrons. The inclusive photon sample is used to select 7 -|- jet events. Jets 
are required to be reconstructed as hadronic tau objects. The true If^T in this sample 
should be zero. The reconstructed If^T corresponds to the instrumental If^T in data. 

The simulation uses Z — > rgT/j process and requires a tight electron and a hadronic 
tau object. The difference between the reconstructed IpT in the simulation minus the 
If^T from neutrinos corresponds to the instrumental Ifix in MC. 

Then the instrumental Ifix is projected to the direction of hadronic tau object, 
shown in Fig. I5.3U1 The distributions in data and MC are different for the longitudinal 
component and the transverse component, respectively. 

To get the uncertainty due to the instrumental IpT-, we "smear" the longitudinal 
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Figure 5.30: Distributions of the instrumental J^t in data using 7+jet sample and 
MC using Z TeTfi sample. Instrumental is projected to the direction of the 
reconstructed leading tau object to get the longitudinal and transverse components. 
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Table 5.13: Number of Z' — > rr events to study the uncertainty in acceptance due to 
the imperfect modeling of the instrumental in MC simulation. The uncertainty is 
obtained from the effect of with/without "smearing" the instrumental in MC to 
that in data. 



component and the transverse component of the instrumental in MC according to 
their differences between data and MC, then add neutrinos back to get the smeared ^t- 
Now we can calculate the uncertainty of the cuts related to J^t by the effect 
with/without smearing the instrumental J^t- Table I^TT^ shows the effect in T^Th chan- 
nel of Z'{m = 300 GeV/c^) sample and Z'{m = 600 GeV/c^) sample. The uncertain- 
ties are (1191-1125)/1125 = 5.9% in Z'(300) sample and (1875-1814)/1814 = 3.4% in 
Z'{600) sample. Taking the larger value, we find that the uncertainty in acceptance 
due to ^ 6%. 
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Chapter 6 
Event Kinematic Selection 

In this chapter we first discuss the trigger paths. Second, we discuss the good run selec- 
tions and the integrated luminosities. Third, in addition to the particle identification, 
we add event kinematic cuts to further suppress backgrounds. Since the kinematic cuts 
need to keep high efficiency for the signals, optimization on the event kinematic cuts 
is necessary. Fourth, there are thresholds in the triggers. The trigger primitives are 
not exactly the same as the offline variables we cut on, and so we need to evaluate the 
marginal trigger efficiencies for selected events. 

6.1 Trigger Path 

For the TgTh selection, we use the "electron plus track" trigger called TAU_ELE. It 
requires an electron in the CEM detector with Et > 8 GeV, px > 8 GeV/c and an 
isolated track with pT > 5 GeV/c. 

For the T^Th selection, there are two "muon plus track" triggers called TAU_CMU 
(TAU_CMX) which requires a CMUP (CMX) muon with pr > 8 GeV/c and an isolated 
track with pT > 5 GeV/c. 

For the ThTh selection, we use the "^t plus tau" trigger called TAU_MET. It re- 
quires LI > 25 GeV and an L3 tau object with Et > 20 GeV, track isolation and 
m(tracks) < 2.0 GeV/cP. 

The TAU_ELE, TAU_CMU, and TAU_CMX triggers are cleaned up by requiring an 
isolated track. The TAU_MET trigger requires only one isolated tau object thus the 
other tau objects in this trigger are not necessarily isolated. 

The track isolation requirement in these triggers is that there is no additional track 
in a 10° to 30° annulus. This track isolation is looser than the offiine tau track isolation 
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with a shrinking inner cone. The detailed descriptions of the tau triggers can be found 
in Ref. 

In addition to selecting the candidate events, there is also an important issue regard- 
ing the jet— > T misidentification background. This fake background is not negligible 
because of the large production rate of jets. Using MC simulation to model all the pro- 
cesses of the fake background is not adequate. We estimate the contribution of these 
events directly from real data. 

For the purpose of estimating jet— > r misidentification background, it is better to 
use those triggers without the isolation requirement in order to have a sample which 
has a larger statistics and is dominated by jet background. 

There is an ELECTR0N_CENTRAL_8 (abbreviated as CELE8) trigger which has 
the same requirement as TAU_ELE but without the track isolation requirement. There 
is also a MU0N_CMUP8 (abbreviated as CMUP8) trigger which has the same require- 
ment as TAU_CMU but without the track isolation requirement. The CELE8 and the 
CMUP8 triggers are dynamically prescaled. A prescale is imposed to reduce the rate 
of a trigger. A fixed prescale under- utilizes the trigger bandwidth when the luminos- 
ity falls during a run. A dynamic prescale is based on the availability of the trigger 
bandwidth, and automatically reduce the prescales as the luminosity falls. 

There is not a corresponding trigger path available for the TAU_CMX trigger. There 
is a prescaled trigger available for the TAU_MET trigger but its prescale is 100 which is 
too big. Thus their jet— > r fake background estimates have to be done with the trigger 
itself. 

6.2 Good Run Selection and Integrated Luminosity 

We use the data samples collected in CDF from March 2002 to September 2003 for 
this analysis. The Good Run List j36j used in this analysis is in the range of the run 
number 141544-168889. 

We use the online initial filtering and the offline periodic classification to decide 
whether a run is good or bad. The former gets rid of obviously bad runs where there are 
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problems with the sub-detectors or the triggers. The latter is based on the classification 
using a large sample in a run, for example, of the J/^' ee, /i/x events which is expected 
to have a very narrow peak, or the photon plus jet events which is expected to have 
very good energy balancing, etc. 

The status of a trigger or a sub-detector is a single bit 1 or 0, which means good 
or bad. The bit 1 or of a trigger is based on whether the deadtime is less than 5% 
and is set by the online run control shift crew. The bit 1 or of a sub-detector at the 
online stage is based on the status of the high voltage, the calibration, the occupancy, 
etc. and is set by the monitoring operator. The bit 1 or of a sub-detector at the 
offline stage is based on, for example, the reconctructed J/^' — > ee, fi/j, mass which can 
tell possible problems in the tracking system, the calorimeters or the muon chambers, 
and is set by the physics groups. 

Here are the details of the requirements on a good run. There are several run 
configurations (trigger tables) when the CDF detector is taking data: test, calibration, 
cosmic, and physics. A good run must be a physics run. At the online stage the 
losses of the beam should be low. The "on-tape" luminosity should be greater than 
10 nb"^. The bits of the LI, L2, L3 trig gers, the calorimeters, the CMU detector, the 
CES detector should be 1. At the offline stage the bits of the calorimeters, the COT 
detector, the CMU and CMP detectors should be 1. The runs after 150145 when the 
CMX trigger updated LI hardware, in addition to the bits above, are required to have 
the online and offline bits of the CMX detector set to 1. 

The total integrated luminosity in the included good runs in the run number range 
141544-168889 is 195 pb"^ However, the good run of the data sample from the 
TAU_CMX trigger starts from 150145 and its integrated luminosity is 179 pb~^; the 
good run number of the data sample from the TAU_MET trigger stops at 156487 and 
its integrated luminosity is 72 pb~^. The TAU_MET trigger was changed after run 
156487 to include L2 two-dimensional track isolation which needs further study. The 
uncertainty in the luminosity measurements is about 6% |27j . 

The integrated luminosity in the data sample from the CELE8 trigger, which is 
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dynamically prescaled, is 46 pb^ . It is calculated by adding the isolated track require- 
ment and comparing its survived number of events with the total number of events in 
the data sample from the TAU_ELE trigger whose luminosity is known. Analogously, 
the integrated luminosity in the data sample from the CMUP8 trigger is found to be 
38 pb-i. 

There were duplicate events incorrectly processed and put in the data samples that 
were later reprocessed. We reprocessed all of the events. To avoid the duplicate events, 
we pick one of them and require that it be a unique event. 

6.3 Selection Criteria 

The event kinematic cuts are designed to further suppress background while to keep 
high signal efficiency. Table 16.11 shows the list of cuts for event selection. We note 
several features of the requirements: 

• The pj, threshold is 25 GeV/c because tau identification is fully efficient at about 
25 GeV/c and it is a high threshold to reduce background. 

• The Ej,, pj,, and pj, thresholds are 10 GeV, 10 GeV/c and 10 GeV/c, respectively. 
(The thresholds in the corresponding triggers are 8 GeV, 8 GeV/c and 8 GeV/c.) 
For T/iT/i, we require the second tau pTp > 10 GeV/c. 

• The If^T cut and the angle cut A(/)(Z — Ipx) < 30" are designed to remove hadronic 
jet backgrounds. They are explained below. 

• We use rriyis > 120 GeV/c^ cut to remove the "irreducible" Z/^* tt back- 
ground. The low mass region with m^is < 120 GeV/c^ is our control region. 

• For the t^t^ selection, we have a cosmic veto j26j . 

• For the ThTh selection, we require the second tau has exactly one track to further 
clean up QCD backgrounds. We will check tau signature by track multilicity on 
the leading tau side. 
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TeTh 




ThTh 


pJp > 25 
E!f, > 10, > 10 
> 15 

A(/){e - ^t) < 30° 
m(e + T + f t) > 120 


plj, > 25 

> 10 

> 15 

A(/)(;U - ^r) < 30° 
m{ii + t + $t)> 120 
cosmic veto 


p^ > 25 
p^? > 10 
> 25 

A(/)(r2 - I^t) < 30° 
m{n + T2 + $t) > 120 
T2 num. track ==1 



Table 6.1: Event kinematic cuts. 



The measured in t^t^ channel needs a muon correction since there is an effect 
of missing energy due to the fact that muons are minimum ionizing partilces. The 
procedure of the muon correction is: first, we subtract the px of a tight muon; second, 
we add muon energy deposits in the calorimeters to avoid counting the same energy 
twice. 

We require > 15 GeV for the TeTh and r^r/j selections. For the ThTh selection, 
we use data from the TAU_MET trigger and we require > 25 GeV to match the 
25 GeV trigger threshold. We could suppress more backgrounds by requiring more 
significant ^t- However, for the signal processes, since there is at least one neutrino at 
each side, there is a chance that the transverse momentum of the neutrinos cancel each 
other, and hence raising thresholds can reduce signal efficiency. We found those 
thresholds are at optimzed points. 

The < 30° cut requires that the significant should follow the e (/x) for the 
TeTh {T^Th) channels and follow the lower pT tau object for the ThTh channel. The 
measured is the vector sum of the neutrinos from the decays of the two taus. Here is 
the example with the TeTh channel which has one neutrino associated with Th and two 
neutrinos associated with Tg. Thus this event topology cut is able to get the most of the 
signals, and to strongly suppress the backgrounds, especially the jet— >^ r misidentified 
fake backgrounds which mostly has a A^ with a random topology. 
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Figure 6.1: LI > 25 GeV trigger efficiency vs. offline for di-tau event. 



6.4 Marginal Efficiency Correction 

We need to include in our estimates of the signal and background rates the effect of 
the triggers. We are concerned, however, only with the effect of the triggers on those 
events passing the offline cuts: the marginal efficiency. 

The TAU_MET trigger for the t^t/j analysis triggers directly on tau object, thus 
there is no marginal trigger efficiency from the TAU side. But there is a marginal 
trigger efficiency from the MET side which is based on a 1 GeV tower threshold for a 
fast calculation at LI while the offiine is based on a 0.1 GeV tower threshold. We 
use the JET20 data sample and mimic the ThTh event topology in the calorimeter by 
requiring one central jet with Et > 25 GeV and at least another one central jet with 
Et > 10 GeV. The LI > 25 GeV trigger efficiency vs. offline for di-tau event 
is shown in Fig. 16.11 The marginal trigger efflciency of the TAU_MET trigger for the 
ThTh analysis is a slow turn-on due to the large trigger tower threshold. 
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The marginal efficiencies of the TAU_ELE and TAU_CMU (TAU_CMX) triggers for 
the TeTh and t^t/i analyses are all at plateau, 

e(TAU_ELE) = 0.92 ± 0.03 (6.1) 
e(TAU_CMU) = 0.85 ± 0.03 (6.2) 
e(TAU_CMX) = 0.92 ± 0.03 (6.3) 

The trigger efficiencies of the electron part, the muon part and the isolated track part 
are calculated by using conversion electrons from 7 — > ee, muons from T/Z — > fifj,, 
and tracks from jet samples, respectively. The details can be found in Ref. [HZI- The 
biggest uncertainty is from the track provided by the XFT trigger, which uses the four 
axial r — (p superlayers (no stereo r — z superlayers) of the COT detector with at least 
10 hits (out of total 12 hits) in each axial superlayer. In the event reconstruction, we 
require at least 3 axial superlayers with at least 7 hits in each axial superlayer, and the 
same configuration for the stereo superlayers. The marginal XFT track finding trigger 
efficiency is found to be a function of pT, the number of prongs, and the different 
run ranges. The overall uncertainty is about 3%. 
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Chapter 7 
Low Mass Control Region 

The low-mass region with rriyis < 120 GeV/c^ is used as the control region to test the 
event cuts and background determination. If we find that the observed and predicted 
event rates agree in the control region, we can proceed to unblind the signal region. 

The main source of events in the control region is from Z/j* tt. The other back- 
grounds include Z/j* ee, Z/^* fifi and jet^ r misidentified fake background. Top 
background tt and di-boson backgrounds such as WW and WZ are negligible because 
their cross sections are two orders of magnitude smaller than Drell-Yan backgrounds 
and their event topology is the opposite of the requirement that a significant fx fol- 
lows the lepton direction. The jet^ r misidentified fake background is not negligible 
because the dijet production cross section is large. 

For the jet— > r misidentified fake background, rather than trying to model all the 
processes that could produce fake events, we estimate the contribution of these events 
from real data which includes any process contributing to the fake background. 

7.1 Drell-Yan Cross Section 

The cross section times branching ratio of the Drell-Yan processes in the mass window 
66 < m < 116 GeV/c^ at ^ = 1.96 TeV is about 250 pb _20_. Fig. O shows the 
mass spectrum and event counts in different mass regions. The Z/j* tt sample has 
377143 events in the mass window 66 < m < 116 GeV/cP which corresponds to a 250 pb 
production cross section. The number of events in a mass window is proportional to 
the cross section in that mass window. For example, the number of events 492000 in 
the mass window m > 30 GeV/c^ gives a cross section 250 x 492000/377143 « 326 pb. 
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Figure 7.1: Drell-Yan mass spectra in different mass regions. 



By tlie same algebra, we get the cross sections in different mass windows: 
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» 326 pb 


a 


B{Zh* - 


-.1+1- 
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(7.1) 
(7.2) 
(7.3) 
(7.4) 



7.2 Drell-Yan Background 

The Drell-Yan backgrounds can be estimated from MC simulation with three pieces: 

Expected MC background = luminosity x (a • B) x acceptance (7-5) 

We just discussed the production cross section, and we have discussed the luminosity 
for each trigger path in Section l6?2l Now we discuss the acceptance and the estimate 
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of the Drell-Yan backgrounds. Table ITTI shows the Drell-Yan background acceptances, 
the apphcation of the trigger efficiencies, the apphcation of the lepton data/MC scale 
factors, and the normalization to the integrated luminosities of the data samples from 
the triggers. 

7.3 Fake Background 

In a "fake" background event a jet is misidentified as a tau. This background is not 
negligible because the dijet production cross section is large. The relative jet— > r 
misidentification rate and the relative tau identification efficiency corresponding to the 
denominator chosen is applied to the denominator tau objects to compute their weight 
for being a jet. We sum up the weights of all the events to get jet— > r misidentified 
fake background estimate in the sample, as described in Section f5. 2. 61 

There is also a probability that, for example, for T^Th channel, a jet is misidentified 
as an electron. But the jet— > e misidentification rate is an order of magnitude smaller 
than the jet— > r misidentification rate. Electron identification requires at most two 
calorimeter towers with EM energy fraction greater than 0.95 and other cuts. Tau 
identification requires at most six calorimeter towers with EM energy fraction less 
than 0.8 corresponding to ^ greater than 0.2 and other cuts. Naively assuming a flat 
distribution between and 6 of the number of towers of jet, and a fiat distribution 
between 0.0 and 1.0 of jet EM energy fraction, we have 

jet^r _ (6 - 0) X (0.8 - 0.0) ^ 

jet ^ e (2 - 0) X (1.0 - 0.95) ^ ' ' 

The electron side is much cleaner than the tau side. It is a good approximation to 
estimate fakes from the tau side. The situation is the same for r^r/i channel. 

There is a subtlety in the fake estimate for ThTh channel. In the data sample from 
the TAU_MET trigger, we order the tau objects in each event by their px- To illustrate 
the subtlety, here we temporately call the leading tau object with the highest pT as ri 
in the case it is a true tau or jeti in the case it is a true jet, and the second tau object 
with a lower px as T2 or jet2. The trigger only requires one isolated tau object. We 
estimate the fake background from the second tau object side which is not necessarily 



96 



source 






— > ee 


1 1 


mass window 


m> 30 


m 


> 30 


m > 30 


a ■ B (pb) 


326 




326 


326 


event 


492000 


398665 


405291 


TeTh (TAU_ELE) 










r(25) + e(10) 


1528 




272 


1 




514 




29 


1 


A6(e - Et) < 30° 


415 




2 





rriina < 120 


405 




1 





trigger efficiency 


373.07 




0.92 


0.00 


lepton scale factors 


351.03 




0.87 


0.00 


normalized (195 pb~^) 


45.36 




0.14 


0.00 


TnTh (TAU CMU) 










r(25) + CMUP it(lO) 


836 







415 


cosmic veto 


836 







415 


$T > 15 


294 







351 


l\(b(u - Et) < 30° 


253 







7 


m,„<j < 120 


248 







4 


trigger cfRcicncy 


226.06 




0.00 


3.65 


lepton scale factors 


191.69 




0.00 


3.09 


normalized (195 pb^) 


24.77 




0.00 


0.48 


TuTft (TAU CMX) 










r(25) + CMX /x(10) 


425 







219 


cosmic veto 


425 







219 


> 15 


150 







181 


Acbia - $t) < 30° 


134 




n 


1 


Tfi.^A^ <C 120 


130 










trigger efficiency 


118.50 




0.00 


0.00 


lepton scale factors 


114.84 




0.00 


0.00 


nornializocl (179 pl^~^) 


13.62 




(J. 00 


0.00 


ThTL (TAU_MET) 










ri(25) + r2(10) 


4264 




1 


9 


f T > 25 


295 










A(f){T2 - ^t) < 30° 


240 










T2 num. track == 1 


185 










rrims < 120 


169 










trigger efficiency 


93.39 




0.00 


0.00 


lepton scale factors 


87.87 




0.00 
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Table 7.1: Drell-Yan background estimates for each channel in the low mass control 
region. 
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isolated. This is able to cover the two cases (a) and (b) out of the total three cases of 
the fake background sources: (a) ri + jet2, (b) jeti + jet2, and (c) jeti + T2. Jeti has 
a lower misidentification rate than jet2 because of its higher pT, so we get c < a and 
a+b ~ a+b+c. The fake estimate from the second tau object side is an approximation. 

The procedure to estimate the jet fake background in the various channels is shown 
in Table UTA We need to define a specific denominator according to data sample from 
the trigger path available, and we need to find out the normalization factors of the 
dynamically prescaled trigger paths. 

The denominators which is up to the electron removal cut ^ > 0.2 and DfrkisoioDeg 
which is up to the 10° track isolation cut are explained in Table I^HI in Section [5.2.31 
Note that the relative jet— > r misidentification rate and the relative tau identification 
efficiency for different denominator samples are different. 

The available dynamically prescaled triggers are discussed in Section lUTTl and their 
integrated luminosities are discussed in Section 16.21 

• The TgT/i channel has a dynamically prescaled data sample from the CELE8 trigger 
path available. There is no trigger cut on the tau objects, so it is ideal for the fake 
background estimate. We apply the cuts up to the electron removal cut ^ > 0.2 
listed in Table 15.11 on the tau objects and use the denominator to estimate 
the fakes. The integrated luminosity of this trigger path is 46 pb~^, thus the 
normalization factor to the integrated luminosity 195 pb^^ of the data sample 
from the TAU_ELE trigger is 195/46 = 4.239. 

• The r^T/i with CMUP muon channel has a dynamically prescaled data sample from 
the CMUP8 trigger path available. There is no trigger cut on the tau objects, 
and we use the denominator to estimate the fakes. The normalization factor 
is 195/38 = 5.132. 

• The T^Tfi with CMX muon channel has to use the TAU_CMX trigger itself for 
the fake background estimate. The tau objects have already been cleaned up 
by the 10° track isolation cut in the trigger. We apply the cuts up to the 10° 
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channel 


TeTh 


CMUP T^Th 


CMX T^Th 




trigger path 
denominator 
norm, factor 


CELE8 
4.239 


CMUP8 
5.132 


TAU.CMX 

DtrklsolODeg 
1 


TAU_MET 
1 


^uj^^* or event 
kinematic cuts 




event 




event 




event 




event 


92.1 
0.903 


2292 
56 


12.4 
0.403 


362 
12 


64.4 
1.649 


379 
30 


106.8 
3.163 


2778 
43 


normahzed 


3.83 ±0.51 


2.07 ±0.60 


1.65 ±0.30 


3.16 ±0.48 



Table 7.2: Fake background estimates in the low mass control region. Uncertainties are 
statistical. 



track isolation cut listed in Table I^TTI on the tau objects and use the denominator 
DtrklsolODeg to estimate the fakes. 

• The ThTh channel has to use the TAU_MET trigger itself. The leading tau object 
is cleaned up by track isolation, but the second tau object is not. We estimate the 
fake background from the second tau object side, and use the denominator D^. 

For each event, we substitute the relative tau identification efficiency and the relative 
jet— > r misidentification rate corresponding to the defined denominator into Eq. (|5.15|) 
if the tau object does not pass the full set of the tau identification cuts, or into Eq. (|5.16|) 
if it does, to calculate the weight to be a jet. 

We sum up the weights of all the events in the sample to estimate the jet background, 
using Eq. ()5.17|) . We then apply the event kinematic cuts and normalize the numbers 
to the luminosities of the data samples of the tau trigger paths. 

The event entries which are integers corresponding to the sum of weights which are 
real numbers are also shown in Table 17.21 The event entries are used to estimate the 
statistical uncertainties. 

There is a systematic uncertainty due to the uncertainty in the jet— > r misiden- 
tification fake rate. The rate used is the average fake rate of the JET samples. We 
use the individual fake rate of the JET20, JET50, JET70, and JETIOO samples to es- 
timate this uncertainty, shown in Table [731 For example, for the r^Th channel, using 
the average fake rate we get an estimate of 3.83; while using the individual fake rates 
from the JET20, JET50, JET70, and JETIOO samples, we get estimates 4.43, 3.25, 
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channel 




CMUP T^Th 


CMX T^Th 




average 


3.83 


2.07 


1.65 


3.16 


JET20 


4.43 


2.23 


1.83 


3.26 


JET50 


3.25 


1.61 


1.35 


2.90 


JET70 


3.03 


1.62 


1.40 


3.01 


JETIOO 


2.94 


1.77 


1.32 


3.20 


syst. err. 


0.89 


0.46 


0.33 


0.26 



Table 7.3: Systematic uncertainties of fake background estimates in the low mass control 
region. 



3.03, and 2.94, respectively. We take the biggest difference, i.e. 3.83 — 2.94 = 0.89 as 
the systematic uncertainty. The fractional systematic uncertainty for this channel is 
0.89/3.83 ~ 20%. The fractional systematic uncertainties of other channels are about 
20% too. 

Combining in quadrature the statistical uncertainties in Table 17^ and the systematic 
uncertainties in Table 1731 we get 

TeTh fake = 3.83 ± 1.03 (7.7) 
CMUP T^Th fake = 2.07 ± 0.76 (7.8) 
CMX T^Th fake = 1.65 ± 0.45 (7.9) 
TftT/, fake = 3.16 ±0.55 (7.10) 

7.4 Cross Check Fake Background 

We perform a cross check on the fake background estimate as follows: relax the tau 
isolation and the lepton isolation, and apply all of the other cuts. The tau isolation 
and the lepton isolation are uncorrelated, thus we can extrapolate from the fake regions 
into the signal region. For example, for rgT/j channel, the signal region A and the 
background regions B, C and D are defined as in Fig. 17.21 and the fake backgrounds in 
A extrapolated = B x D/C. 

Unfortunately, we can only cross check the fake background for TeT^ channel using 
the data sample from the CELE8 trigger path and possibly for r^r/i with CMUP muon 
channel using the data sample from the CMUP8 trigger path. Neither sample has 
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Figure 7.2: Using the uncorrelated tau isolation and electron isolation to estimate fake 
background for TgTh channel. 

isolation in the trigger. There is no such sample for t^t/j with CMX muon channel. 
There is a sample without isolation for t/jT/j channel, but its prescale is 100 which is 
too large for this exercise. 

Due to the statistics in the region B, C and D, this cross check can only be done 
for the TeTh channel in the low mass control region. The numbers in region B, C and 
D are 12, 142 and 13, respectively. When we extrapolate to region A, we find that 

A = B X D/C = 12 X 13/142 = 1.099 (7.11) 

The normalization factor is 4.239, thus we get TeT^ fake extroplated = 1.099 x 4.239 = 
4.66. This is in good agreement with 3.83 it 1.03 obtained by summing up the weights 
of tau object being a jet. It does give us confidence in the method of the jet— > r 
misidentified fake background estimate. 



101 



7.5 Uncertainties in Control Region 

The statistical uncertainty and the systematic uncertainty of Dreh-Yan background 
estimate include 

• statistical uncertainty, 

• a ■ B uncertainty, 2%, aside from luminosity uncertainty (see Ref. j20j). 

• trigger efficiencies (see Section . 

• lepton scale factors (see Section [5.3.31 for r scale factor, Section [5.4.21 for e scale 
factor, and Section [5.5.21 for /i scale factors), 

• uncertainty, 6% (see Section [^?B|) . and 

• luminosity, 6% (see Ref. [23 )• 

The statistical uncertainty and systematic uncertainty of the jet— > r misidentified 
fake background estimate are discussed in Section [7.31 

We combine the t^t^ CMUP muon channel with a luminosity 195 pb~^ and the 
T^Tfi CMX muon channel with a luminosity 179 pb^^ into one channel, simply called 
the T^T/i channel. The observed events in TeTh, r^Th and ThTh channels are 46, 36 and 
8, respectively. 

Table [7711 shows the summary of the control sample in low mass region for 195 pb~^ 
(72 pb~^ for ThTh). The total background estimate is 99.27 it 12.55, dominated by the 
source from Z/^* tt as expected. The observed number of events, 90, in the control 
region is in good agreement with this prediction. Fig. 17. 31 - 17. 5[ show the distributions 
of each channel in the low mass control region. The observed distributions in the data 
are in good agreement with the predicted distributions. 
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Source 








Total 


Z/7* ^ TT 


45.36 ± 6.84 


38.39 ± 5.72 


4.19 ±0.77 


87.94 ± 12.38 


Z/"f* — > ee 


0.14 ±0.14 








0.14 ±0.14 


Z/7* ^ /x^ 





0.48 ±0.25 





0.48 ± 0.25 


Jet— > r 


3.83 ± 1.03 


3.72 ±0.88 


3.16 ±0.55 


10.71 ± 1.46 


Expected 


49.32 ± 6.94 


42.59 ± 5.85 


7.35 ± 0.95 


99.27 ± 12.55 


Observed 


46 


36 


8 


90 



Table 7.4: Number of expected events for each channel and each source, compared with 
the number observed, in the control region m^j^ < 120 GeV/c^. 
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Figure 7.3: Distributions of the T^Tfi channel in the control region for data (points) and 
predicted backgrounds (histograms). 
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Figure 7.4: Distributions of the r^r/i channel in the control region for data (points) and 
predicted backgrounds (histograms). 
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Figure 7.5: Distributions of the r^r/i channel in the control region for data (points) and 
predicted backgrounds (histograms). 
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Chapter 8 
High Mass Signal Region 

The high mass region with rriyis > 120 GeV/c^ is the signal region. First we calcu- 
late signal acceptance, then we estimate the backgrounds. The main backgrounds are 
Z/7* TT, Z/j* — > ee, Z/j* fifi which can be estimated from MC simulation, and 
the jet— > r misidentified fake background which can be estimated from data, as in the 
control region. 

8.1 Signal Acceptance 

Table ISTTl shows the procedure to measure the signal acceptances in each channel for the 
new vector particle decaying to two taus, using Z' tt events. For example, for the 
TeT/i channel, we match the offline tau object and electron object with the and Te by 
requiring the separation angle be less than 0.2 radian, apply the event kinematic cuts, 
multiply the number of accepted events by the trigger efficiency and the lepton scale 
factors, and calculate the overall acceptance. Since the mass of the Z' is unknown, we 
calculate the signal acceptance as a function of its mass. Only five mass points (120, 
180, 300, 450, 600) GeV/c^ out of total twelve mass points (120, 140, 160, 180, 200, 250, 
300, 350, 400, 450, 500, 600) GeV/c^ are shown in Table EH The signal acceptances of 
the T^Th channel with a CMUP muon and of the t^t/j with a CMX muon are combined 
into one signal acceptance for the t^t/j channel. The total acceptance is a combination 
of the acceptance of the rgT/j, the t^t/^, and the t/jT/j channels. The signal acceptances 
are shown in in Fig. 18.11 

Table 18.21 shows the procedure to measure the signal acceptances in each channel 
for the new scalar particle decaying to two taus, using A — > tt. We set tan /3 = 20 as 
a representative value of tan /?. Similarly, since the mass of A is unknown, we calculate 



106 



the signal acceptances as a function of mass, as shown in Fig. 18.21 

8.2 Drell-Yan Background 

The largest portion of the production cross section for the Drell-Yan backgrounds is at 
the Z boson resonance peak, about 91 GeV/c^. However the events in the high mass 
signal region are mostly from the high mass Drell-Yan tail. To model the high mass tail 
better, we need more statistics in MC simulation at that region. To achieve this, we 
break the generation level mass into two exclusively separated regions: 30 < m < 100 
GeV/c^ and m > 100 GeV/c^, and simulate them separately. The production cross 
sections in these two regions are about 315 pb and 11 pb, respectively (see Section [7?T|) . 
Therefore we have a low-mass sample and a high-mass sample for each Z/j* 
source. 

Table 18.31 shows the procedure to estimate Drell-Yan backgrounds. We apply the 
event kinematic cuts on the MC samples, multiply the number of surviving events by the 
trigger efficiencies and the lepton scale factors, normalize to the integrated luminosity 
195 pb~i (179 pb-i for the TAU.CMX trigger, 72 pb^^ for the TAU_MET trigger), 
and combine the estimate for the low-mass Drell-Yan sample and the estimate for the 
high-mass Drell-Yan sample. 

8.3 Fake Background 

The procedure to estimate the jet— > r fake background is similar to what we have 
done for low mass control region in Section 17.31 The trigger path, the luminosity 
normalization factor, the denominator tau object definition, and the sum of the weights 
of tau objects being a jet in the high mass signal region are exactly the same as those in 
the low mass control region. The only one difference is this cut: rriyis < 120 GeV/c^ for 
the low mass control region, while rrims > 120 GeV/c^ for the high mass signal region. 
Now we repeat the same procedure, as shown in Table 15^ The event entries which are 
integers corresponding to the sum of weights which are real numbers are also shown. 
The event entries are used to estimate the statistical uncertainties. 
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576.1 
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1015.4 
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861.1 


acceptance (%) 
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0.682 
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2449 


3415 


3932 
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> 25 


277 


940 


2037 


2888 


3383 


AoiCro — Et) < 30° 


242 


832 


1679 


2244 


2459 


T2 num. track == 1 


185 


653 


1335 


1789 


2043 


rrins > 120 


31 


526 


1282 


1768 


2028 


trigger efficiency 


21.2 


388.3 


1023.7 


1469.1 


1716.9 


lepton scale factors 


19.9 


365.3 


963.2 


1382.3 


1615.4 


acceptance (%) 


0.020 


0.365 


0.963 


1.382 


1.615 


channels combined 












acceptance (%) 


0.045 


0.907 


2.673 


3.756 


4.377 



Table 8.1: New vector particle Z' — rr signal acceptance, for each channel, as a 
function of the Z' mass. 
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Tabic 8.2: New scalar particle tt signal acceptance, for each channel, as a function 
of the A mass. 
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Figure 8.1: Signal acceptance of a new vector particle Z' — > rr in each channel, as a 
function of the Z' mass. 




Figure 8.2: Signal acceptance of a new scalar particle A 
function of the A mass. 
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Table 8.3: Drell-Yan background estimates for each channel in the high mass signal 
region. 



Ill 



channel 


TeTh 


CMUP T^Th 


CMX T^Th 


Th-Th 


trigger path 
denominator 
norm, factor 


CELE8 
D_xi 
4.239 


CMUP8 
D_xi 
5.132 


TAU.CMX 
D.trklsolODeg 
1 


TAU_MET 
D_xi 
1 


^uj^^* or event 
kinematic cuts 




event 




event 




event 




event 


92.1 
0.068 


2292 
13 


12.4 
0.006 


362 
1 


64.4 
0.152 


379 

4 


106.8 
0.282 


2778 
12 


normahzed 


0.29 ±0.08 


0.03 ± 0.03 


0.15 ±0.08 


0.28 ±0.08 



Table 8.4: Fake background estimates in the signal region. Uncertainties are statistical. 



channel 


TeTh 


CMUP T^Th 


CMX Tf^Th 


ThTh 


average 


0.29 


0.03 


0.15 


0.28 


JET20 


0.18 


0.03 


0.16 


0.31 


JET50 


0.23 


0.04 


0.15 


0.23 


JET70 


0.31 


0.03 


0.14 


0.25 


JETIOO 


0.28 


0.03 


0.13 


0.22 


syst. err 


0.11 


0.01 


0.02 


0.06 



Table 8.5: Systematic uncertainties of fake background estimates in the signal region. 



There is a systematic uncertainty due to the uncertainty in the jet^ r fake rate. 
The rate used is the average fake rate of the JET samples. We use the individual fake 
rate of the JET20, JET50, JET70, and JETIOO samples to estimate this uncertainty, 
as shown in Table 

Combining in quadrature the statistical uncertainties in Table lSTH and the systematic 
uncertainties in Table we get 

TeT,, fake = 0.29 ±0.14 (8.1) 
CMUP T^Th fake = 0.03 ± 0.03 (8.2) 
CMX r^r,, fake = 0.15 ±0.08 (8.3) 
ThThi&ke = 0.28 ±0.10 (8.4) 

8.4 Uncertainties in Signal Region 



We summarize all of the systematic uncertainties in the high mass signal region in this 
section. Some of these are due to statistical uncertainties on the various backgrounds 
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Z' TT 


m = 120 


m = 180 


m = 300 


m = 400 


m = 600 


MRST72 / CTEQ5L 


1.047 


1.029 


1.021 


1.006 


1.002 


MRST75 / MRST72 


0.951 


0.980 


0.983 


0.995 


0.993 


CTEQ6L1 / CTEQ6L 


1.006 


1.006 


1.003 


0.999 


1.002 


CTEQ6M / CTEQ5L 


1.035 


1.023 


1.021 


1.008 


1.004 


PDF uncertainty 


7.7% 


4.2% 


3.4% 


1.1% 


0.8% 



Table 8.6: PDF uncertainty. 



due to limited Monte Carlo or other statistics. Others come from separate external 
studies as indicated. And in this section, we combine the r^Th with CMUP muon 
channel and the T^Th with CMX muon channel into one single t^t/j channel. 

The systematic uncertainty in the Drell-Yan and new particle signal rates due to 
the imperfect knowledge of the parton density functions (PDF's) is calculated 
by comparing the acceptance change ratio for various PDF's. The CTEQ5L is used in 
PYTHIA. We add in quadrature the difference between MRST72 to CTEQ5L, MRST75 
to MRST72, CTEQ6L1 to CTEQ6L, and CTEQ6M to CTEQ5L PDF's. The MRST72 
and MRST75 compare the effect of varying as on the PDF. The CTEQ5L set is leading 
order, and the CTEQ6M sets are next to leading order but at the same value of a^. 
Using Z' TT, this is shown in Table We take 8% as a conservative number. 

We are careful to identify the correlated and the uncorrelated systematic uncertain- 
ties. The correlated uncertainties include the uncertainties of the PDF, the integrated 
luminosity, the e, /x, r scale factors, the ^t, and the jet^ r fake rate. Table IHTTI lists 
the uncertainties, their magnitude, and the affected channels. (When uncertainties are 
correlated we assume a 100% correlation.) 

The Z' — > TT and A ^ tt signal acceptances and the systematic uncertainties are 
listed in Table 18.81 - 18.91 The acceptance itself reflects the effects of trigger efficiency 
and the lepton scale factors. The uncertainties include the contributions from 

• statistical uncertainty (MC statistics), 

• PDF uncertainty (this Section), 

• trigger efficiencies (see Section . 
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Uncertainty 


Magnitude (%) 


Affected Channels 


PDF 


8 


all 


integrated luminosity 


6 


all 


e scale factor 


4 




// scale factor 


5.5 




r scale factor 


10 


all 




6 


all 


jet— > r fake rate 


20 


all 



Table 8.7: Systematic uncertainties, in percent, and the affected channels. 



m(Z') 


TeTh (%) 


r^r^ (%) 


ThTh (%) 


combined (%) 


120 


0.012 ±0 


004 


0.013 ± 0.004 


0.020 ± 0.005 


0.045 ± 0.009 


140 


0.084 ± 


015 


0.088 ±0.015 


0.105 ±0.020 


0.278 ± 0.043 


160 


0.213 ±0 


035 


0.151 ±0.025 


0.206 ± 0.038 


0.571 ±0.086 


180 


0.308 ±0 


049 


0.234 ±0.037 


0.365 ± 0.066 


0.907 ±0.136 


200 


0.453 ±0 


070 


0.351 ±0.054 


0.476 ± 0.085 


1.280 ±0.190 


250 


0.727 ±0 


111 


0.548 ± 0.083 


0.776 ±0.137 


2.052 ± 0.303 


300 


0.975 ± 


148 


0.735 ±0.110 


0.963 ±0.170 


2.673 ± 0.394 


350 


1.098 ±0 


167 


0.826 ±0.124 


1.144 ±0.202 


3.068 ± 0.452 


400 


1.239 ±0 


188 


0.966 ±0.144 


1.308 ± 0.230 


3.512 ±0.517 


450 


1.395 ±0 


211 


0.979 ±0.146 


1.382 ± 0.243 


3.756 ± 0.553 


500 


1.537 ±0 


232 


1.148 ±0.172 


1.431 ±0.252 


4.116 ±0.604 


600 


1.572 ±0 


237 


1.190 ±0.178 


1.615 ±0.284 


4.377 ± 0.644 



Table 8.8: Uncertainties of // — > — > rr signal acceptance (SM coupling). 



• lepton scale factors (see Section [5.3.31 for r scale factor. Section 15.4.21 for e scale 
factor, and Section 15.5.21 for scale factors), and 

• ]pT uncertainty (see Section EH) . 

The systematic uncertainties on the Drell-Yan backgrounds and the jet— > r misiden- 
tified fake backgrounds are listed in Table 18.101 The systematic uncertainties on the 
Drell-Yan backgrounds incorporate the effects of 

• statistical uncertainty (MC statistics), 

• PDF uncertainty (this Section), 

• a ■ B uncertainty, 2%, aside from luminosity uncertainty (see Ref. [20]). 
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m{A) 


TeTh (%) 


T^Tf, (%) 


ThTh (%) 


combined (%) 


120 


0.020 ± 0.005 





018 ± 0.005 


0.029 ± 0.007 


0.067 ±0.012 


140 


0.113 ±0.019 





082 ± 0.015 


0.126 ±0.024 


0.321 ±0.050 


160 


0.284 ± 0.045 





213 ± 0.034 


0.324 ± 0.058 


0.822 ±0.123 


180 


0.474 ± 0.074 





376 ± 0.058 


0.445 ± 0.079 


1.295 ±0.191 


200 


0.603 ± 0.093 





485 ± 0.074 


0.660 ±0.117 


1.748 ±0.259 


250 


0.889 ±0.135 





703 ±0.106 


0.972 ±0.172 


2.564 ±0.379 


300 


1.174 ±0.178 





846 ±0.127 


1.132 ±0.199 


3.151 ±0.463 


350 


1.254 ±0.190 


1 


004 ±0.150 


1.356 ± 0.239 


3.614 ±0.532 


400 


1.411 ±0.213 


1 


101 ±0.165 


1.485 ± 0.261 


3.996 ± 0.588 


450 


1.460 ± 0.220 


1 


135 ±0.170 


1.573 ± 0.277 


4.168 ±0.614 


500 


1.649 ± 0.249 


1 


177 ±0.176 


1.561 ±0.275 


4.386 ± 0.644 


600 


1.758 ±0.265 


1 


300 ±0.194 


1.684 ±0.296 


4.742 ± 0.696 



Table 8.9: Uncertainties oi gg ^ A ^ tt signal acceptance (tan/? = 20). 



Source 





rr 


z/i*- 


ee 


z/r- 


MA* 


Jet- 


> T 



TeTh 



0.56 ±0.11 
0.16 ±0.14 


0.29 ±0.14 



T^,Th 



0.50 ±0.10 


0.50 ±0.26 
0.18 ±0.09 



Expected 1.01 ±0.24 1.18 ±0.30 0.64 ±0.13 2.83 ± 0.46 



ThTh 



0.36 ± 0.08 



0.28 ±0.10 



Total 



1.42 ± 0.23 
0.16 ±0.14 
0.50 ±0.26 
0.75 ±0.19 



Table 8.10: Uncertainties of backgrounds in signal region, 195 pb (72 pb for ThTh). 



trigger efficiencies (see Section IHT!)) . 



• lepton scale factors (see Section [5.3.31 for r scale factor, Section 15.4.21 for e scale 
factor, and Section 15.5.21 for fj. scale factors), 

• uncertainty (see Section [^^ . and 

• luminosity, 6% (see Ref. [17]). 

The systematic uncertainties on the jet^ r misidentified fake background incorporates 
the effects of 

• statistical uncertainty (see Section . and 

• systematic uncertainty due to jet— > r misidentification rate (see Section . 
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Chapter 9 
Results 

9.1 Observed Events 

After unblinding the signal region, we observe four events in TeT^ channel, zero events in 
T^T/i channel, and zero events in r^r/i channel. The numbers of background events esti- 
mated and observed are in Table IHTTl Fig. 19. II shows the rrims distribution. Fig. 19.21 - 19.51 
shows the event displays of the four events observed in TeTh channel. 

9.2 Experimental Limits 

Since we observe no excess, we proceed to calculate the 95% confidence level (CL) 
upper limit on the cross section times branching ratio for new particle production using 
a Bayesian procedure described in Ref. [HS] . 

We need to combine multiple search channels and incorporate both uncorrelated 
and correlated systematic uncertainties. For each channel i, the integrated luminosity, 
the signal acceptance, the expected background events, and the observed events are 
denoted as Lj, e^, 6j, and rii, respectively; the uncorrelated uncertainties of the signal 



Source 


TeTh 


T,,Th 


ThTh 


Total 


Z/7* TT 


0.56 ±0.11 


0.50 ±0.10 


0.36 ±0.08 


1.42 ±0.23 


Z/7* ee 


0.16 ±0.14 








0.16 ±0.14 


z/r ^ 





0.50 ±0.26 





0.50 ± 0.26 


Jet^ T 


0.29 ±0.14 


0.18 ±0.09 


0.28 ±0.10 


0.75 ±0.19 


Expected 


1.01 ±0.24 


1.18 ±0.30 


0.64 ±0.13 


2.83 ± 0.46 


Observed 


4 








4 



Table 9.1: Number of expected events for each channel and each source, and number 
of observed events, in the signal region. 
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Figure 9.1: Distribution of visible mass {rrims) for data (points) and predicted back- 
grounds (histograms) in the signal and control regions. The upper plot is in linear scale. 
The lower plot is in log scale. 




Figure 9.2: Event display TeTh candidate run=152669 event=629080 m^,js=148 GeV/c^. 




Figure 9.3: Event display rgT/j candidate run=153693 event=815662 mms=^'^^ GeV/c^. 
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Figure 9 



.5: Event display rgT/j candidate run=162252 event=612118 myis=^'^^ GeV/c^. 
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acceptance and the expected background events are denoted as f^i and fbi, respectively. 
The correlated uncertainties of the integrated luminosity, the signal acceptance, and the 
expected background events are denoted as gi, Qe, and g^, respectively. (Note that the 
/ factors carry i indices and the g factors do not.) With a signal cross section asig-, the 
expected number of event ^Uj in each channel can be written as 

^Jii = {l+ gL)Liasig{l + fei){l + ge)ei + (1 + /m)(1 + 9b)bi (9.1) 

where the / and g factors are in a form 1 + x thus relative systematic uncertainties. 
We define a likelihood which is the product of the Poisson probabilities of observing Ui 
events in each channel, 

C{n\as^gAe)=Y[/:{n,\^l,) = l[^^^-^ (9.2) 

i j 

where the overbars indicate that the variables are arrays carrying an i index. We use 
a Monte Carlo method to convolute the effects of the systematic uncertainties using 
Gaussian prior probability density functions for the / and g factors. For an evaluating 
point of the asig, we sample the / and g factors within their Gaussian widths around 
a central value of zero, calculate the /ij and the £(nj|/Xj) for each channel, and average 
the resulting likelihood C{n\asig,b,e). Using Bayes' Theorem, we then construct a 
probability density function for the signal cross section. 



C{n\asig,b,e)P{asigj 
/o°° ^(.nla'^tg^ b, e)Pia'^^g)da'^ig 



n'^s^gin, b, e) = ^rrrr^^ r . , (9-3) 



with a prior probability density function P[asig) which expresses the subjective "degree 
of belief" for the value of the signal cross section. The 95% CL upper limit (T95 is 

obtained by solving this integral equation 

0-95 

V{asig\nAe)d(Jsig = Q.'ib (9.4) 

We assume a uniform prior in the signal cross section up to some high cutoff; the value 
of the cutoff has no significant influence on the 95% CL upper limit. 

We thereby extract the experimental 95% CL upper limit of a ■ B for models using 
vector boson and scalar boson, respectively. The results are listed in Table 19.21 and 
shown in Fig. 19.61 These are the generic limits for gg ^scalar ~^ '^'^ ^^'^ ~^ 
"'^vector ~^ '^'^ which can be interpreted in various models. 
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mass 


vector 


scalar 


(GeV/c2) 


limit (pb) 


limit (pb) 


120 


122.294 


87.338 


140 


18.884 


17.899 


160 


9.446 


6.996 


180 


6.066 


4.229 


200 


4.185 


3.187 


250 


2.637 


2.192 


300 


1.999 


1.764 


350 


1.757 


1.540 


400 


1.537 


1.396 


450 


1.441 


1.330 


500 


1.296 


1.290 


600 


1.237 


1.174 



Table 9.2: The 95% CL upper limits on vector and scalar particle production and decay 
to tau pairs. 
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nix [GeV/c^] 

Figure 9.6: Upper limits at 95% CL on a{pp X)B{X tt) for vector and scalar 
boson, as a function of mass. 
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Figure 9.7: Upper limits at 95% CL and theoretical predictions of a{pp Z')B{Z' 
tt). The excluded region is the region with m{Z') < 394 GeV/c^. 



9.3 Exclusion Regions 

Now we can put the theoretical predictions on high mass tau pair production discussed 
in Section f2.3l and the experimental 95% CL upper limits together. We take the region 
where the theoretical prediction is bigger than the upper limit to be excluded at 95% CL. 

For reference, this analysis would thus exclude at 95% CL a Z' with standard model 
couplings having a mass of less than 394 GeV/c^, as shown in Fig. MJ\ For the MSSM 
pseudoscalar Higgs boson A, this analysis is not sensitive to exclude a region yet. 
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Chapter 10 
Conclusions 



We have performed a blind search for high mass tau pairs using data corresponding 
to 195 pb^^ of integrated luminosity from Run II of the Tevatron, using the CDF 
detector. In the high-mass region with rrims > 120 GeV/c^, we expect 2.8 it 0.5 events 
from known background sources, and observe 4 events in the data sample. Thus no 
significant excess is observed, and we use the result to set upper limits on the cross 
section times branching ratio to tau pairs of scalar and vector particles as a function of 
mass, shown in Table 19.21 and ploted in Fig. 19.61 
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Appendix A 



The Structure of the Standard Model 



The fundamental constituents of matter in Nature are fermions: leptons and quarks, 
with interactions specified by the gauge symmetries SU(3)cxSU(2)ixU(l)y in the 
framework of the Standard Model (SM). 

Why are the fermions in an electroweak doublet? Why are left-handed fermions in 
a doublet, and right-handed fermions in a singlet? What tells us that quarks have color 
degrees-of-freedom? Why must quark doublets be paired with lepton doublets? Here 
we come to a brief review of how the structure of the SM emerged. A good introduction 
can be found in Ref. j39j . 

The relationships among the fermions are interpreted from the interactions they 
experience, namely the cross sections and decay widths measured, calculated, and mea- 
sured ... an interplay of experimental inputs and theoretical constraints. The objective 
is to unify the different interactions. 

Charged Current 

Let us recall Fermi's theory jlUj of charged current (CC) weak interaction for four 
fermions, e.g. the crossed /3-decay, ep rei^e, 



P 




n 



e 



The amplitude (matrix element) of this process can be written as 



(A.l) 
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where Gp is, Fermi's constant and the charged currents for the fermion fields are 



(A.2) 



The next advance came after the discovery that CC violates parity maximally j41| 
and the V-A theory of the weak interaction [IJj was proposed. Only left-handed 
fermions, which are projected by a V-A operator ^(1 — 75), appears in CC. 

helicity allowed 




M 



V 


- 1/2 


V 


V 


-1-1/2 


V 


V 


-1-1/2 


X 


V 


- 1/2 


X 


[ = 


V2 ^ 





,1 



1 



(A.3) 
(A.4) 



After the introduction of quarks 1^ for understanding the classification of the 
hadrons, it was natural to re-write the hadronic part of CC in terms of quark fields. 
The transition u — > d occurs via CC, with the other two quarks in the nucleon being 
spectators. 




(A.5) 



There was an inconsistency found in the value of the Fermi constant as de- 
termined from /3-deay and the purely leptonic muon decay. This lead Cabibbo to the 
hypothesis that the quark states in CC are not the physical states (eigenstates of mass), 
but rather a quantum superposition of the physical states. 



d 



weak 



cos 9 c sin Oc 
^ — sin 9c cos 9c 



( 



d 



(A.6) 



mass 
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where 6c is Cabibbo angle, thus the Fermi constant is replaced hy Gp cos 9c ■ This idea 
was generalized to the case of three quark generations in terms of the CKM (Cabbibo- 
Kobayashi-Maskawa) matrix 



\ / weak 



V 



Vub 

Vcd yes Vcb 
Vtd Vts Vtb 



(A.7) 



mass 



Glashow proposed the intermediate vector boson model (IVB) in 1961 [1^ and the 
form has been incorportated in the SM. The basic idea is to replace the four fermion 
interaction by the exchange of a massive charged boson VF^, e.g. v^e~ —>■ lyefJ-' , (a) 
four fermion interaction, (b) the IVB model: 




The matrix element can be written as 




Fermi 



^Gf jCc 
V2 ^ 



l^jccy 



(A. 



IVB 



9 jCC 
V2 



m 



w , 



(A.9) 



Comparing Eq. HA.8|) with Eq. (|A.9|) . substituting g = e/ smOw and a = e^/(47r) 



and using experimental values: a = 1/137, Gp = 1.166 x 10 GeV , sin 9 



w 



0.22 



(sin^ 9w was first determined from the NC/CC cross section ratio in neutrino scattering 
where NC is the neutral current interaction explained below) , this leads to the prediction 
for the W mass: 



mw 



37.3 



SGp I sin^vF 
This may be compared with the experimental value 0: 



79.5 GeV/c^ 



(A.IO) 



mw = 80.425 ± 0.038 GeV/c^ 



(A.ll) 
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The interdediate bosons, along with the bosons explained below, were discovered 
at CERN in 1983 US). 



A Doublet in Weak Isospin Space 

We write the left-handed leptons in a weak isospin SU(2)i doublet and the right- 
handed leptons in a singlet, for example, 

/ \ 

L= \,eR (A.12) 

V J 

The generators of the SU(2)i transformations are = ^r*, where r* are Pauli matrices. 
The charge raising opertator r+, the charge lowering operator r~, and the original 
are 



-{ri+iT2) 



(n - iT2) 

1 
-1 




(A.13) 
(A.14) 
(A15) 



The currents can be written as 



J, 



Jf, 



T^elfi-{1 - 75)e = l'el,ieL = Lj^T+L 

1 _ _ - 1 



(A.16) 
(A.17) 
(A.18) 



These can be combined into an isospin triplet of currents 



(A.19) 



The weak isospin invariance implies that the SU(2)i invariant Lagrangian to describe 
the interaction between the W bosons and the current J with a coupling g is of the 
form 

C = g3^'. (A.20) 
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Hence a neutral IVB should exist, coupling to J^. Since the electromagnetic current 
J'^ is parity conserving, 

= e{eR-f^eR + eLl^^ei) (A.21) 

Whereas has a V-A structure. cannot be directly identified with the electromag- 
netic current, nor with the photon. 



Neutral Current 

Next came the inputs from the neutral current (NC) interactions. NC were discov- 
ered by the Gargamelle Collaboration at CERN in 1973 |46j . v^q v^q. 
CC . NC 




■■'^ >-hadrons 
The matrix element can be written as 




? >-hadrons 

P- 



^ ^ 8^ l^jNCy jNC 



with NC in the form 



+ E 

/ 



1 



(A.22) 



(A.23) 



l = e,fi,T; f = l,q; q = u,d,s,c,b,t 

The neutrino part has a V-A structure. The lepton/quark part has parity violation 
{C{ ^ 0), but not maximally {C{ ^ Cl). Universality of NC and CC requires p = 1, 
later predicted in the SM. 

We can write the NC interactions in terms of IVB, e.g. f^e~ — > f^e~, (a) four 
fermion interaction, (b) the IVB model: 








propagator J 


\z 
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The matrix element can be written as 

Mfy^B - ^-^r i^] (j'"'')' (A-25) 

COS Ow \ fn^ ) cos Qw ^ ' 

Comparing Eq. (|A.24|) with Eq. ()A.25|) . and assuming universahty of the charged and 
neutral currents (/J = 1), this gives the prediction for Z mass: 

i^9^\''' 1 raw 79.5 p 2 

"^^= r^T^ ^ Z— = ^ Z— = ^ = 90GeV/c A.26 

y 8Gf y yTJcost^vy ^cosdw cost^iy 

This may be compared with the experimental value 0: 

mz = 91.1876 ± 0.0021 GeV/c^ (A.27) 



Flavor Changing Neutral Current 

The flavor changing neutral current (FCNC) interaction is strongly suppressed 
(a) CC, (b) FCNC: 

^- ^ 

_[7t+ 




s ^ w-i ^ e 




J^C' = _^ ^'t^' = uu + dd cos^ + ss sin^ 60 + (sd + ds) sin 6*0 cos (A.28) 

FCNC 

BR(i^+ ^ 7r+7r°) = (21.13 ± 0.14)% (A.29) 

BR(K+ ^ vr+e+e") = (2.88 ± 0.13) x 10"^ (A.30) 

The GIM (Glashow-Iliopoulos-Maiani) mechanism [17j proposed that quarks must be 
paired in doublets. This naturally solved FCNC. In addition, c quark was predicted 
and later discovered |48j . 

J^*^ = UU + CC + d'd' + ss' = uu + CC + dd + ss (A. 32) 
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A Triplet in Quark Color Space 

The quarks in the spin-| baryons are in a symmetrical state of space, spin and flavor 
degrees of freedom, e.g. 

A^^ = uuu, Q~ = sss (A. 33) 

However the requirements of Fermi-Dirac statistics imply the total antisymmetry of the 
wave function. The solution was the introduction of the color degree of freedom, with 
indices as red (r), green (g), and blue (6). 



(A.34) 



V '^^ / 

One of the tests of the number of charged fundamental constituents is provided by 

R = f ^ (A.35) 

The virtual photon emitted by the e^e~ annihilation will excite all kinematically ac- 
cessible qq pairs from the vacuum. 



At low energy where only the u, d and s quarks are available, in the absense of color 
degree of freedom, we expect 

If quarks have three colors, 

R = Z{el + el + el) = 2 (A.38) 

For energies above 10 GeV, c and b quarks are available, 

R = 2,{el + el + el + el + eg) = y (A.39) 

The color triplet model is excellently supported by data, see the "cj and R in e^e~ 
Collisions" plots in the Section "Plots of cross sections and related quantities (Rev.)" 
in PDG 0. 
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Pair Quarks with Leptons 

Some classical symmetries, known as anomalous symmetries H^l are broken by 
quantum effects. The requirement for an anomaly- free theory ]50' is that: 



where the sum is over all quarks and leptons. For example consider the two doublets, 



Cancellation of anomalies requires that quark doublets must be paired with lepton 
doublets. The SM identifies a generation in a natural way by identifying the doublet 
containing the heaviest charged lepton with the doublet containing the heaviest quarks 
(and so on), but one could in principle associate any quark doublet with any lepton 
doublet and call that a generation, because there are no interactions between quarks 
and leptons in the SM. What needs to be guaranteed is that the number of quark and 
lepton generations must be equal. 




(A.40) 






(A.41) 
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Appendix B 

Gauge Symmetry & Spontaneous Symmetry Breaking 

The interactions between the fermions and the vector bosons in the Standard Model 
(SM) are uniquely specified by requiring the theory, i.e. the SM Lagrangian, invariant 
under gauge transformations which are local and involve transformations varying from 
point to point. Some of the standard texts are listed in Ref. jSlj . 

A symmetry indicates a deeper relationship among the elementary particles with 
a further unification of the interactions and makes the form of a Lagrangian more 
compact. Symmetry dictates design and plays the central role in the direction to find 
the simplest model. 

Gauge Symmetry 

Let us take electromagnetism as an example and consider the Lagrangian for a free 
fermion field ^{x). 



This is invariant under a global U(l) phase transformation which is space-time inde- 
pendent and is illustrated in the left plot in Fig. IB. 11 



where Q is the charge or the U(l) quantum number of the fermion. For example, the 
charge assignment for u quark, d quark, i/g, and e are +2/3, -1/3, 0, and -1, respectively. 

We are going to construct an invariant Lagrangian under a local, i.e., gauge, U(l) 
phase transformation which is space-time dependent and is illustrated in the right plot 
in Fig. EH 



Co = ^(2;)(i7^9^ - m)'^{x) 



(B.l) 



(B.2) 



(B.3) 



134 







/« 






■ \ 


t 












/* 








/ 


V 




\ 










/* 






t 


I 












/* 


1 t 






\ 


\ 












1 ^ 




\ 


I 


4t^ 



Global Local 



Figure B.l: Global and local transformations. 

The partial derivative 9^ in Eq. ()B.1|) spoils the invariance. We need to form a gauge- 
covariant derivative which will have the simple transformation property, 

D^,^{x) ^ e-^«^(^)L»/,^(x) (B.4) 

so that the combination ^Dp^ is gauge invariant. To achieve this, we enlarge the 
Lagrangian with a new vector gauge field and form the covariant form as 

= {d^,^ieQAp)^ (B.5) 

where e is a free parameter which eventually will be identified as the coupling of the 
gauge field to the fermion field. The transformation property in Eq. ()B.4|) will be 
satisfied if the gauge field A^ix) has the transformation property 

Ap{x) ^ Alpix) = Ap{x) + ^d,,9{x) (B.6) 

Note that the coupling of the gauge field (photon) to any fermion field is determined 
by its transformation property under the symmetry group. This is usually referred to 
as universality. Also note that photon is massless because an Af^A^^ term is not gauge 
invariant under this transformation. 

To make the photon field a truely dynamical variable we need to add a kinetic term 
to the Lagrangian involving its derivatives. The simplest gauge-invariant term with a 
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conventional normalization is 




(B.7) 



where 





(B.8) 



Terms with higher powers are omitted in order that the theory be renormalizable. We 
notice that photon does not have self-coupling because it does not carry a charge. 
Now we arrive at the gauge-invariant QED Lagrangian 



Most remarkably, if one demands the symmetry be local, one is forced to include the 
electromagnetic field, and hence, light. Recall that there are four Maxwell equations. 
While here we just require "gauge symmetry" and electromagnetism is determined. 
This illustrates how physics becomes simpler. 

Non-Abelian Gauge Symmetry 

Yang and Mills extended the gauge principle to non-Abelian symmetry j52j . Con- 
sider the simplest case isospin SU(2). Let the fermion field be an isospin doublet. 



(B.9) 




(B.IO) 



The free Lagragian 



(B.ll) 



is invariant under the global SU(2) transformation 



-iT-e 



(B.12) 



where 9 = {91,62,0^) are the SU(2) transformation parameters and T = ^ are the 
SU(2) generators with r = {ti,T2,t^) the Pauli matrices satisfying 



[Ti, Tj] = ieijkTk k = 1,2,3 



(B.13) 



with eijk the structure constants for SU(2). 
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It is easy to check that two successive SU(2) transformations do not commute be- 
cause the generators do not commute and this is why SU(2) is called a non-Abelian 
symmetry, in contrast to an Abelian symmetry such as U(l) where two successive U(l) 
transformations commute. 

Under the local symmetry transformation 

^ ^'{x) = e-*'^-^(^)*(x) (B.14) 

the partial derivative in Eq. ()B.11|) spoils the invariance. To construct a gauge- 
invariant Lagrangian we follow a procedure similar to that of the Abelian case: 

• We form a gauge-covariant derivative 

D^^ix) ^ e-'^-^^'-'^D^^ix) (B.15) 

by introducing vector gauge fields j4^, z = 1,2,3 (one for each group generator) 
and a coupling g 

D^^ = {d^ + igT-A^)^ (B.16) 
and defining the transformation property for the vector gauge fields as, 

Al^4=Al- e'^'e^Af, + -^d,e' (B.17) 

The gauge fields are massless because an A'^^A^^ term is not gauge invariant, 
similar to an Abelian field. But, the second term is clearly the transformation for 
a triplet representation under SU(2), thus the fields carry charges. 

• Then we add a gauge invariant kinetic term for the gauge fields 

£a = --/l^F"^^'' (B.18) 

where 

= d^Al - d,A'^ - ge'^'^Aj.At (B.19) 

The third term shows that the gauge fields have self-coupling because they carry 
charge, in contrast to an Abelian field. 
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We arrive at the complete gauge-invariant Lagrangian which describes the interac- 
tion between the gauge fields and the SU(2) doublet fields, 

C = *Z7^(5^ + igT ■ A^)* - m** - ^F'^^F''"' (B.20) 

Generalization of the Yang-Mills theory to a higher group SU(N) with A'" > 3 is 
straightforward . 

SU(3)cXSU(2)LXU(l)y 

The structure of the gauge symmetries in the SM is SU(3)cxSU(2)iXU(l)y. For 
a particular fermion its quantum field is a product of factors, 

^ f space-time \ f spin \ f U(l)y \ / SU(2)^ \ f SU(3)^ 

I factor I \ factor I \ factor / y factor I \ factor 

(B.21) 

Each factor has some labels, coordinates, or indices. The orthonormality of the quantum 
field holds separately for each factor. Since the gauge bosons of one of the symmetry 
groups do not transform under the other gauge symmetries in the product of groups, the 
gauge invariant Lagrangian may be simply written as a sum of the terms of individual 
groups. The gauge symmetric Lagrangian in the framework of Yang-Mills theory is 

jC-symmetric = ^j'' (^8,, + igijB^ + ig2Tmf, + igsX^C;^ ^ (B.22) 

where the eight and A", the three WJ^ and T*, the one and Y arc the gauge 

bosons and generators corresponding to the SU(3)c color, the SU(2)l weak isospin, 
and the U(l)y hypercharge gauge symmetries, respectively; gi are the gauge couplings; 
and 

= dfj^Bi, — diiB^j, (B.23) 

= d^Wl-d,Wl + g2e''^Wl,W^ (B.24) 

= d^Gl-d^Gl + g^r'-GlGl (B.25) 

with e^^^ and /"'"^ the structure constants for SU(2) and SU(3). 
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At this stage, all of the gauge bosons and fermions are massless. The explicit mass 
terms break gauge invariance. For gauge bosons, the exptected mass terms 

m'^wW^W (B.26) 

plus similar terms for the others, arc clearly not invariant under gauge transformations 
WJ^ — s- WJ^ = WJ^ — e^^^O^Wj^ + ^dij,9^. This is true for any gauge theory. For fermions, 
using the left- and right-handed projection operator Pl and Pr, the mass term can be 
written as 

= vti^PlPl'^ + rn^PnPR^ 

= mi^R^L + ^L'^R) (B.27) 

In the SM, left-handed fermions are in SU(2) doublets and the right-handed fermions 
are in SU(2) singlets, thus they transform dilfcrcntly. The '^r'^'l and 'if^'^R terms are 
not SU(2) singlets and would not give an SU(2) invariant Lagrangian. 

However the description that all of the gauge bosons and fermions are massless is 
not true in Nature. We need to 

(a) generate the masses of the leptons and quarks; 

(b) generate the masses of the W'^, W~ , and Z'^ weak vector bosons; 

(c) but also keep the photon and gluon massless. 

In other words, the SU(3)c will be kept precise, and the gluon will remain massless. 
We need to break SU(2)iXU(l)y down to U(1)bm, resulting in mixing between the 
and fields, and non-zero masses for three of the gauge bosons {W^ and The 
photon (A) remain massless, due to a residual U(1)em gauge symmetry that remains 
unbroken. 

Spontaneous Symmetry Breaking 

The solution in the SM is to add a spontaneous symmetry breaking (SSB) term 
into the symmetric Lagrangian "by hand" . The Lagrangian will remain symmetric but 
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the physical vacuum does not respect the symmetry. In this case, the symmetry of the 
Lagrangian is said to be spontaneously broken. 

= ^symmetric + ^SSB (B.28) 

The assumption to construct CssB is that the universe is filled with a scalar field, 
called Higgs field. One real scalar field could solve (a). One complex field could solve 
(a) and create one massive vector boson. To achieve (a), (b) and (c), the minimum 
requirement of the Higgs field is two complex fields arranged in a doublet in the SU(2) 
space and carries U(l) hypercharge +1 (electric charge Q = + ^ is +1 and for the 
upper and lower component, respectively), but is a singlet in color space. 



/ 



(B.29) 



_ 1 91+ 192 

Under a SU(2)2,xU(l)y gauge transformation, the doublet transforms as 

The scalar field can be given gauge invariant terms in CssB- the kinetic term 
required by gauge invariance, the Higgs potential including a mass- like term and a 
self-interaction term, and the Yukawa coupling between the doublet and a particular 
fermion ^. 

CsSB = {D^(t>)\D^'^) - y ((/>) - Cvukan^a (B.31) 



with 



= d^ + mlB^ + ig2Tmi (B.32) 



= /x2</.t(/, + A ((/>V) (B.33) 

jO-Yukawa = 5/^0^ (B.34) 

Spontaneous symmetry breaking of the Higgs potential [1] is possible by assuming 
< (also a positive A to possess a stable vacuum). This is shown in Fig. IB.21 The 
minimum of the Higgs potential shifts (in field space) from (p = to 

= ^ + '/'i + 0i + = ^ = (B.35) 
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Figure B.2: Spontaneous symmetry breaking of Higgs potential. 



The field thus acquires a non-zero vacuum expectation value (VEV). Choosing {(p^) = v, 
we expand about v, 

01 + i<t>2 
V + H + i(f)i 

with (Ir' = v + H. Any SU(2) doublet can be written as 




(B.36) 







(7 X 



(B.37) 



By applying the gauge symmetry oi CssB under the transformation of the Higgs doublet 
in Eq. (|B.30() . the algebra can be simplied by "gauging away" three of the four real 
degrees of freedom of the Higgs doublet with (/)^, c/)^, c/)^ = 0, 

/ 



1 

71 



(B.38) 



V 





V + H{x) 

This is called the unitary gauge. On the other hand, the physical quantities are inde- 
pendent of the choice of gauge. This indicates these degrees of freedom are unphysical. 



141 



Gauge Boson Mass 

The generators of the SU(2)j;, transformations are 
matrices. 

/ 

, I.I 




^r*, where r* are Pauh 



(B.39) 
(B.40) 
(B.41) 



We write exphcitly 



1 1 1 wj^ 

g, B^ + g2Tmi=gi-B^ + g2- ^ 

2 2 ^\wl + tWl 



wl-iwl 



(B.42) 



We then substitute Eq. HB.38|) and Eq. (|B.42|) into the kinetic term and the Higgs 



potential of CssB in Eq. ([B.^ljl . After some algebra the tree- level mass terms for the 
H field and the gauge bosons are present. The unphysical scalars reappear as the 
longitudinal polarizations of the weak bosons. 



{D^<t>)\D^^ct>)-V{<t,) 



4m 



(B.43) 
(B.44) 



+ (2' {92Wl-giB^)[g2W^>'-giB^') (B.45) 



and many other interaction terms. The fields are defined as the electric charge 

eigenstates. The SSB has mixed the and gauge bosons with the weak mixing 
angle Ow- 

1 „-oX 



cos Ow — sin 9w 
sin Ow cos Ow 



(B.47) 
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tan 6w = — 

92 



Now we can read out the tree- level masses, 



(B.48) 



(B.49) 
(B.50) 
(B.51) 
(B.52) 

Using mw = i^^Y^'^ in Eq. (fOn|) with Fermi's constant Gp = 1.166 x 10^^ GeV"^ 



rriH 


= vV2X 


mw 


" 2 


mz 


92 

= V 

2 cos 6w 




= 



we can estimate the VEV of the Higgs field: 



The quantity 



v = ^ = {V2Gf)-'/' « 246 GeV 

92 



mw 



(B.53) 



(B.54) 



mz cos 6w 

is the universality parameter of the neutral current interactions and the charged current 
interactions. It is predicted to be one at tree level in the SM, thus provides a test of 
the SM realization of SSB compared to other models. Any deviation from p = 1 would 
be an important signal of new physics. 



Eletroweak Unification 



Substituting the physical state of in Eq. I|B.46|) and Z^, in Eq. (|B.47|1 
into the electroweak interaction in the covariant derivative term in Eq. HB.22|) . and 
using Y = 2{Q — T^), we can identify the weak CC, weak NC, and electromagnetic 
interactions. 

Y 



-^7^ 91 {Q- T^) + <72 [T'WI, + T'^W^ + T^W^ 
+ r+W-)l^ (weakCC) 



£1. 

V2 



92 



COS 



-sin^OwQ) ^ (weakNC) 



(electromagnetic) 



(B.55) 
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Comparing the electromagnetic part with the — term of Cqed in Eq. HB.9|) . 
this imphes the unification relation: 

e = (72 sin 9^ = gi cos O^j (B.56) 



Yukawa Coupling 

Now we check the fermion masses. The structure of the lepton fields, for example, 
of the first generation is 

L = 1^ j , CR (B.57) 

The Higgs field is an SU(2) doublet. This makes it possible to write an SU(2)-invariant 
interaction of the fermions with the Higgs field, i.e., the Yukawa coupling term in 
Eq. (|R34l) . which can be written as 

C-Yukawa = geL(f)eR + geeRcj)^ L (B.58) 

= Qei X X) \\^]+9e[x](x X 



Here Lcj) is an SU(2) invariant. Multiplying by the cr does not change the SU(2) 
invariance. The second term is the Hermitian conjugate of the first. The coupling 
is arbitrary because it is not specified by the gauge symmetry principle of the theory. 
After SSB by substituting (f) with Eq. (|B.38|) . and using clgr + crcl = ee, we get 

C-Yukawa = ^ (eLCR + eRei) + ^ {e-LCR + eRCi) H 

TTl 

= rrieee H -eeH (B.59) 

V 

We have identified the fermion mass as mg = Thus the theory can now accom- 
modate a non-zero fermion mass. The second term says that there is a lepton-Higgs 
coupling We notice that there is no mass term occured for neutrinos, = 0. By 
assumption the theory contains no right-handed neutrino state i^r, therefore a term 
analogous to Eq. (|B.58|) cannot be written that will lead to a mass term vrvl. And 
this implies neutrinos do not interact with H. 
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The structure of the quark fields, for example, of the first generation is 

/ 



Q 



L 



, UR, dR (B.60) 



UL 

Since the structure of the right-handed quark is different from the lepton case, there is 
a subtlety in writing down the Yukawa coupling term. We know cj) is an SU(2) doublet, 
then so is 

4) = iT^d)* = ' I (B.61) 





This is true for any SU(2) doublet. Since (p has hypercharge Y = +1, (p has Y = — 1, 
and for each state, Q = T"^ + Y/2 is still satisfied. After SSB, cj) becomes 

(B.62) 

The SU(2)-invariant Yukawa coupling for the quarks can be written as 

C-Yukawa = 9dQL4>dR + guQL4>UR + h.C. (B.63) 

After SSB by substituting cp with Eq. HB.38() . (p with Eq. HB.62|1 . and using QlQr+QrQl = 
qq, we get 

^Yukawa = "^^^ + + "71 

= rtiuuu + rtiddd -\ -uuH -\ -ddH (B.64) 

V V 

Again the quark masses can be accommodated, but are arbitrary parameters. They 
have to be provided by experiment. The last two terms describe the interaction of u 
and d quarks with H. 

The procedure can be copied for the second and third generations with e ^ fi,T 
and with u ^ c,t and (i — > s,b. Since H interacts with a coupling proportional to m/, 
it couples most strongly to the heaviest generation. 

CKM Matrix 

The spaces we have been working on are an internal quantum phase space called 
gauge space of fermions, and an internal field space of the Higgs potential. The logic 
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line is gauge symmetry + SSB. Let us write down the SM Lagrangian HB.28|1 explicitly 
by combining Eq. (|B.22() and Eq. ()B.31() . This time is not for a particular fermion 
only. There are three generations of fermions in the SM. We will sum up all of them. 
Once we do that, there is a new internal space: generation space. The eigenstates of 
the fermions in gauge space could be not the eigenstates of the fermions in generation 
space which are the physical mass eigenstates we observe in experiment. 



^ — ^symmetric ~l~ ^SSB 



(B.65) 



+ 



d^, + igi\B^ + ig2Tmi, 



2 



, covariant ) 
{,Lsymm, Gk) 

kinetic) 
{I' SSB, V{4>)) 

-Hfgf^ f(t>^ f {LsSB, Yukawa ) 

We collect all of the terms for the fermions after SSB: the kinetic and QCD terms in 
Eq. (|B.22|) . the mass and Higgs coupling terms in Eq. HB.59|) for the leptons and in 
Eq. ()B.64|) for the quarks, and the weak CC, weak NC and electromagnetic terms in 
Eq. ()B.55|) . We simplify the notation for a fermion field j as /. The part of the SM 
Lagrangian for fermions is given by 



Cf= Eff{ip-mf-^H)f 
-eEfQlfl'^fA^ 



(Higgs) 
(QCD) 
(QED) 



(B.66) 



- ^ E/ /V {T' - sin2 ewQ) fZ, (weak NC) 

- ^ E/ /V(T+VF+ + T~W-)f (weak CC) 
We denote the gauge eigenstate triplets in the generation space as 



UL 



/ \ 

^J'L 

V / 
/ \ 

UL 
CL 



/ \ 

en 

V ^« / 
/ \ 

UR 
CR 



(B.67) 
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1 


dL 


\ 




/ 




\ 


dL = 




SL 








SR 






\ 


bi 


) 




I 


bR 


J 



and denote the rotations from the gauge eigenstates to the mass eigenstates as unitary 
matrices Lg, Re, Lu, Ru-, Ld, and Rd such that 

RefiR 

RuUR (B.68) 
RddR 

Because ah of the neutrinos in the SM are massless, they are degenerate in the mass 
eigenstates, namely we cannot teU the differences among the mass eigenstates. We set 
the rotation for neutrinos as a unit matrix denoted as I^- 

First we check the QED part in Eq. (fR66|) to see if there is any change under the 
rotations, 

= -ej:Qfh'fA, 

f 

= - eQf{hl''h + hYfR)Af„ withf = e,u,d 

^ -eQf[hl^LlLfiL + h7''RlRfiR)A^ (B.69) 

= - eQf {hr^L + hl^fR) 

where we have let pass 7^ forward in Eq. (|B.69() because the former rotates in 

the generation space and the latter is in the spinor space. Since the unitary rotation 
matrices give L^jLf = I and R^jRf = I, the electromagnetic interaction is diagonized 
in both the gauge eigenstates and the mass eigenstates. 

The same result holds for the Higgs, QCD, and weak NC parts in Eq. (|B.66|) for 
the same reason. For the weak NC, this is called the GIM mechanism |47| . The flavor 
changing neutral currents (FCNC), e.g. s ^ d decay "off-diagonal" in the generation 
space, are strongly suppressed. On the other hand, the FCNC rare decays are very 
interesting because they are possible probes for new interactions. 



Ul LuUL, \1R 

dL LddL, dR 
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Now we check the weak CC in Eq. HB.66|) . For leptons, we have 



CC 



92 

V2 



uj^blW + h.c. 



(i>4Le7^eLTy^- + h.c. 
^(ur^LW- + h.c. 



(B.70) 



where we have let pass backward in ()B.70|) . With I^Lf. acting backward on the 
vector of the degenerated neutrino mass eigenstates, we just go back to the original form, 
and the leptonic weak CC interactions are diagonized in both kinds of the eigenstates. 

So far, the distinction between the gauge eigenstates and the mass eigenstates has 
been seen to have no apparent effect. However, mixing between generations does man- 
ifest itself in the system of the weak CC for quarks. By convention, the quark mixing 
is assigned to the down-type quarks, 



-CC 



^ (nL^dLW; + h.c) 
(iiLrVdLW- + h.c. 



(B.71) 



92 ( 



where 



(B.72) 



Thus the down- type quark gauge states participating in the transitions of the weak CC 
are linear combinations of their mass eigenstates. For three generations, it is called the 
CKM (Cabibbo-Kobayashi-Maskawa) matrix 5 . The SM does not predict the content 
of V . Rather its matrix elements must be extracted from experiment. 

/ 



V / weak 



V 



(B.73) 



Vud Vus Vub 
Vcd Vcs Vcb 
Vtd Vts Vtb 

Any 3x3 complex matrix has 18 paramters. The quark mixing matrix V ^ being 
the product of two unitary matrices, is itself unitary, V'^V = 1, and this eliminates 



V/ 



mass 
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9 paramters. The rest of 9 parameters can be identified with 3 rotation angle, and 6 
phase angles with 5 of them eliminated by rephasing the relative quark phase angles 
in Eq. HB.73|) and leaving 1 global phase angle. So the actual total number of free 
parameters is 18 — 9 — 5 = 4, which includes 3 rotation angle and 1 phase angle. 
The "standard" parametrization of the CKM matrix advocated in PDG |3j is 



V 



\ 



C12C13 

-S12C23 - Cl2S23'Sl3e 



S12C13 



-i(5i3 



'013 



C12C23 - Sl2S23'Sl3e 



(B.74) 



S23C13 

S12S23 - Cl2C23Sl3e'"" -C12S23 - Sl2C23'Sl3e''"" C23C13 

In this equation, dj = cos Oij and Sij = sin Oij , with i and j labeling the genera- 
tions. The interpretation is that if 6ij vanishes, so does the mixing between those two 
generations. For example, in the limit 823 = ^13 = 0, the third generation decouples 
and it reduces to two generations with 612 identified as the Cabibbo angle. 

The complex parameter in phase angle goes into the weak charged interaction terms 
u'y^^PlVdWf^, and from quantum theory we know that the Hamiltonian will not be 
invariant under time reversal, or equivalently, CP. So this induces CP violation. 

The magnitude of the complex matrix element in the CKM matrix presently mea- 
sured is 



^ 0.9739 - 0.9751 0.221 - 0.227 0.0029 - 0.0045 ^ 



V 



(B.75) 



0.221 - 0.227 0.9730 - 0.9744 0.039 - 0.044 
0.0048 - 0.0140 0.037 - 0.043 0.9990 - 0.9992 
Here we discuss some of the immediate consequences. For top quark, with Vtb ~ 0.999, 
we have 

BR(t ^ Wb) ^ 100% (B.76) 

For bottom quark, with Vet ~ 0.04 ten times larger than ~ 0.004, it mostly decays 
by 6 ^ Wc. Then W can decay to ev, fiv, tv, ud, and cs with a color factor 3 for each 
quark decaying mode. The width of b decays is Pf, (9yJ,G|.m^)/(1927r^). This gives 



rv ^ 5^ ' ~ 



0.4 



(B.77) 



So the life time of b is about two and half times longer than the life time of r because its 
decay can only happen by the rotation from the mass eigenstates to the weak eigenstates 
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and the magnitudes of the matrix elements for this rotation are small. 



Couplings to Fermions 



For convenience- 



we repeat Eq. ()B.66|) here. 



dp 



(T^-sin^^H^Q)/^, 



(Higgs) 



(weak CC) 



(weak NC) 



(QED) 



(QCD) 



(B.78) 



We can read out the couplings to the fermions in the SM as follows: 



• The Higgs coupling for — > // is — . 



• The QCD coupling for g ^ qq is 



^A". (For the electroweak interactions of the 



quarks ^/Z/W/H — > ijcqc, the effect of the color charge is that the probabilities, 
i.e., the decay widths are multiplied by a constant color factor Nc = 3, rather 
than that the couplings appearing in the amplitudes are multiplied by the color 
generator A"^. This is because that ^/Z/W/H are colorless and the number of 
color combinations of gc<?c is fixed to be three.) 

• The electromagnetic coupling for 7 ^ // is eQf. 

• The neutral weak coupling for Z^ — > // is cof Ivf (-^'^ ~ ^^^'^ ^wQf) for left-handed 
fermions and ( —22RjJ}^—2^\ for right-handed fermions. 

\ cos Ow J 

• The charged weak coupling is and this only applies to left-handed fermions. 
We notice that the coupling for —>■ Ivi is while the coupling for —>■ qq' 
should be multiplied by a quark mixing element in the CKM matrix and it becomes 




These results are summarized in Table ITU in Section [2.11 
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Appendix C 
How to Calculate Cross Section 



We are concerned about the resonance production of tau pairs in the SM i.e. pp 
7*/Z TT. This is a good example to see how event generator |5] |53j works by using 
Monte Carlo simulation. 

At pp collider, the production of any process starts from parton interaction. A 
proton is made of quarks and gluons and can be written as 

proton = uud + uu + dd + --- + g + g + --- (CI) 

'^■^ V ' ^ ' 

valence sea gluons 

The probability density for a given parton i in a proton carrying momentum fraction 
X and being "seen" in an interaction by an intermediate boson with energy scale Q is 
characterized by a function fi{x,Q), called the Parton Density Function (PDF) |22j . 
The momentum density of a parton is its PDF multiplied by its momentum fraction 
and is expressed as xfi{x,Q). An example of parametrization is shown in Fig. IC.ll 
The differential cross section for 12 34 can be written as 

= ^'^02 ^slj^ ,^ ts'oF '^^(Pl +P2 -P3 -P4)rfXidx2/l(3:i)/2(x2) I-^I?2^34 

is (27r) 2^3 (27r) 2^4 

(C.2) 

where s is the parton center-of-mass energy squared, pi (Ei) is the momentum (energy) 
of the ith particle, xi^2 are the fractions of the momenta of the incoming beam particles 
carried by the incoming partons, fi{xi) are the PDF's with an implicit dependence on 
the energy scale of the interaction, and Yl,spins l-^li2^34 matrix element squared 

for the process averaged over the spins and colors of the incoming particles and summed 
over the spins and colors of the outgoing particles. 

First we consider the phase space. We perform the integral over the three-momentum 
of and reexpress the integral over the momentum of p^ in terms of the magnitude of 
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X 

Figure C.l: Proton's parton density functions. 

the three-momentum p in the parton center-of-mass frame and the angle with respect to 
the beam 6 and the azimuthal angle (f). Then we make a transformation using s = xiX2S 
with s the pp center-of-mass energy squared and we get dx2 = ds/{sxi). After some 
algebra, the differential cross section becomes 

d(7 = ^^-^^d cos ed(l)ds^xif i{xi)x2f2ix2) Yl l-^li2^34 (^-3) 

^ '^^ spins 

The angular part can be uniformly generated with < (/) < 2it and — 1 < cos 9 < 1. 
The momentum fraction dxi/xi part can be transformed to Inxi and then uniformly 
generated. For the distribution over s, we impose a minimum value of s. There are two 
types of distributions to be smoothed in order to converge faster for the Monte Carlo 
simulation. One type is a power law distribution with a > 1 which is the rise 
in the cross section due to the photon exchange at small center-of-mass energies. The 
other type is the Breit-Wigner resonance due to the Z boson exchange with a mass m 
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and a width T, 

Is/s p = lnxi ^ J dp uniformly 

§ p = s(i-") ^ [dp uniformly (C.4) 

/L„ (s-m^f+m^r^ P^tan-i(l^) Jdp uniformly 

Second we consider the matrix element which is the interesting part. With a non- 
constant matrix element, the distribution is expected to deviate from the pure phase 
space distribution. Further, compared with the distributions described by the SM, there 
are probably deviations in the distributions in real data due to some unknown matrix 
elements of new physics. The effects shown in cross section could be an enhancement 
or a suppression, a new resonance, changes in the angular distributions, a divergence 
or a cancellation by interference, etc. A good deal of particle physics consists of the 
measurements and the interpretations of such effects in cross section. For the SM 
process — > 7*/Z ^ rr, we have 

El-^l?2^34= {i-m^,){i-ml){\g^^\^ + \g^^\^) (C.5) 

+2m3mine{g^^ g^^* + g^^g^^*} 

where t = {pi — Ps)^, u = {pi — ^4)^, 771,3^4 are the masses of the outgoing tau particles. 
In the center-of-mass frame using = ^ [s — {ni^ + m^)"^] [s — {m^ — m^)^] , the value of 
i can be expressed as t = mg — s^^'^^E^ —Pcm cos 9), and the value of u can be expressed 
as n = ml — s^/^(£'4 — pcm cos 6). The couplings are defined to be 

ab \ ^ 9in9out /p p\ 

where the sum runs over 7*/Z the intermediate gauge bosons with mass Q/mz and 
width O/r^, and g^^^ is the coupling of the gauge boson to the incoming partons and 

L Ft 

g^^i is the coupling of the gauge boson to the outgoing tau particles. The couplings to 
the fermions in the SM are listed in Table ITU 

To summarize, the major parts for generating an event include: (a) generating 
randomly the incoming partons and incorporating the PDF's, (b) generating randomly 
the kinematic variables which describe the event in the phase space of the final particles, 
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and (c) calculating the matrix element. Now we can put the parts together and get the 
weight for an event by multiplying all of the factors. The weight is in GeV^^ and we 
need to convert to picobarn with a conversion constant 3.89379 x 10^ GeV^ pb. 

After generating a large sample of events, we can fill the weights of the events into 
a histogram, for example, a one-dimensional histogram of s which is the invariant mass 
of the tau pairs. The differential cross section versus the invariant mass of the tau pairs 
can be obtained by dividing the histogram by the number of events generated and the 
size of the bins. The result is shown in Fig. 12.41 in Section [2. 31 
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Appendix D 
Separation Angle under Boost 



The calculable case is to boost the simplest phase space, i.e. a two-body decay, from 
the rest frame to the lab frame, as shown in Fig. ID. II The two final particles are back- 
to-back in the rest frame. The separation angle a of the two final particles in the lab 
frame can be parametrized as a function of 0* the polar angle in the rest frame, which 
has an equal probability to be any value between 0° and 90°, and the boost 7. 

Let us consider two massless final particles, e.g. two photons from a vr'' decay with 
a mass m, an energy E, and a boost 7 = E/m. We boost the four-momentum of pi 
from the rest frame to the lab frame. 
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We denote the angle between pi in the lab frame and the direction of the boost as 61. 
We have 

sin 6** , 
sin6'i = — ^ = (D.2) 

\l sin^ e* + (7 cos e* + ^Ji^ - 1)2 
We denote the angle between p2 in the lab frame and the direction of the boost as 02- 
By substituting Q* with 6** -|- vr, we have 
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/sin^ Q* (-7 cose* - 1)^ 

Now we can calculate the separation angle ce. 
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For 9* not too small and 7 ^ 1, we get an approximation for small a, 

12 , , 

a « — — X - (D.5 

sm 0* 7 

For fixed 9* (not too small) values, the functions are shown in Fig. ID. 21 For large 
boosts, the smearing by 9* is small, thus the correlation between the separation angle 
and the boost (energy) is very strong. 

Since 9* has an equal probability to be any value between 0° and 90°, the probability 
that the separation angle stays between the curve for 9* = 30° and the curve for 9* = 90° 
is three times larger than the probability that the separation angle stays between the 
curve for 9* = 10*^ and the curve for 9* = 30°. The effect is very obvious. We use 
Monte Carlo simulation to check the same plot, as shown in Fig. ID. 31 It confirms that 
the simplest case of two-body decay is indeed calculable and the correlation between 
the separation angle and the boost (energy) is very strong. 

For the more complicated phase spaces such as those of tau's hadronic decays, the 
calculation is very hard. But Eq. ()D.5|) is still a good hint. We need to use Monte Carlo 
simulation to get the distribution, which is shown in Fig. 15.41 in Section [5.2.21 
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Boost Y (= E/m) 



Figure D.2: Separation angle vs. boost, calculated in 6* slices. 
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Figure D.3: Separation angle vs. boost, Monte Carlo distribution. 
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